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SECTION  1 
INTRODUCTION 


1.1  BACKGROUND 

"k 

The  objective  of  the  proposed  PLACES  barium  release  experiment  is 
to  help  determine  the  effects  of  propagation  disturbances  produced  by  striated 
plasma  regions  on  the  performance  of  phase-coherent  satellite  systems,  with 
emphasis  on  the  Global  Positioning  System  (GPS).  To  relate  experimental 
results  as  directly  as  possible  to  GPS  receiver  performance,  one  would  prefer 
tests  using  actual  GPS  L-band  navigation  signals  if  interesting  levels  of 
propagation  disturbances  could  b^  produced  in  the  field  experiment.  In  this 
context  "interesting"  means  disturbance  levels  comparable  to  those  expected 
in  nuclear  environments.  Previous  calculations  of  propagation  disturbances 
(Reference  l-l)  predict  that  GPS  L-band  signals  are  likely  to  be  driven  into 
saturated  Rayleigh-like  scintillation  in  striated  nuclear  environments.  It 
is  interesting  to  note  that  a  Wideband  Satellite  pass  at  Ancon  on  30  March 
1977  shows  that  saturated  L-band  scintillation  is  occasionally  observed  in 
the  equatorial  ionosphere  (Reference  1-2).  Significant  UHF  scintillation  was 
produced  by  STRESS  barium  releases  (Reference  1-3).  However,  it  is  un¬ 
likely  that  significant  scintillation  will  occur  at  L-band  in  barium 
clouds.  Thus  it  is  natural  to  consider  the  use  of  somewhat  lower  fre¬ 
quencies  in  the  PLACES  experiment. 


*  PLACES  is  the  acronym  for  Position  Location  and  Communication  Effects 
Simulat inn. 
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Because  of  the  wide  bandwidth  GPS  pseudo-noise  spread-spectruir. 
navigation  signals,  selection  of  an  appropriate  frequency  requires  considera¬ 
tion  of  frequency-selective  propagation  mechanisms  including  angular  scatter¬ 
ing,  dispersive  effects,  fading  spectra  and  scintillation  intensity.  These 
disturbances  are  functions  of  the  spectra  of  the  ionization  irregularity 
structures,  as  well  as  propagation  geometry,  frequency,  and  signal  spectrum. 
For  example,  integrated  phase  variance  is  sensitive  to  the  effective  outer 
scale  size,  whereas  angulc.  scattering  is  sensitive  to  the  effective  inner 
scale  size  (which  may  be  limited  by  the  receiver  antenna  aperture).  Scintil¬ 
lation  intensity  is  sensitive  to  the  magnitude  of  the  irregularity  spectrum 
near  the  Fresnel  size,  which  depends  on  frequency  and  propagation  geometry. 

The  point  here  is  that  the  potential  propagation  disturbances  are  sufficiently 
numerous  and  their  relationsliips  to  plasma  properties  sufficiently  complicated 
to  warrant  detailed  calculations  to  establish  proper  scaling  with  carrier 
frequency  and  signal  spectrum. 
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1.2 


APPROAC!! 


The  goal  of  this  report  is  to  provide  pre- experimental  predictions 
of  the  GPS  X-set  receiver  performance  after  signal  propagation  through  a 
disturbed  barium  environment.  This  goal  is  accomplished  here  through  the 
following  steps: 

1.  A  detailed  statistical  model  for  a  striated  barium  cloud  is 
developed  based  on  back-propagated  results  provided  by  ESL, 

Inc.,  from  measurements  made  during  the  STRESS  test  series. 

As  discussed  in  Section  2,  measurements  made  during  three 
aircraft  passes  beneath  the  striated  STRESS  barium  clouds  are 
analyzed  to  obtain  the  integrated  phase  power  spectral  density 
(PSD)  including  rms  phase  fluctuations,  outer  scale  size,  and 
spectral  index. 

2.  This  statistical  model  is  utilized  to  obtain  frequency-selective 
multiple  phase  screen  (MPS)  calculations  of  signal  propagation 
over  signal  bandwidths  I'equired  for  the  GPS  pseudo-noise 
spread- spectrum  navigation  signals.  Carrier  frequencies  range 
from  150  to  250  Mllz. 

3.  The  previously  existing  GPS  X-set  receiver  simulation  was 
modified  to  include  the  received  signal  spectrum  and  the 
receiver  code  correlation  operation.  As  discussed  in  Section 
3,  these  modifications  are  the  only  changes  in  the  receiver 
simulation  required  for  frequency-selective  effects. 

4.  The  frequency-selective  propagation  calculations  called  out 
in  Item  2  above  are  then  used  as  direct  input  to  the  X-set 
receiver  simulation  to  measure  the  receiver  performance  after 
signal  propagation  through  a  barium  environment.  In  Section  4 
the  <^requency-selective  propagation  results  and  the  X-set 
receiver  performance  characteristics  are  discussed. 


It  is  important  to  note  that  all  results  presented  in  this  report 
are  based  on  the  detailed  X-set  receiver  simulation  discussed  in  Reference  1-1. 
This  simulation  is  a  direct,  sampled-data  digital  implementation  of  the  GPS 
X-set  receiver  initially  designed  by  Cahn  (Reference  1-4)  and  described  in 
Reference  1-5.  With  the  exception  of  a  reduced  carrier  frequency  and  restora¬ 
tion  of  the  250/sec  AFC/Costas  iteration  rate,  it  is  assumed  that  there  are  no 
changes  in  the  X-set  design.  Thus  even  with  a  carrier  frequency  as  low  as 
150  MHz,  it  is  assumed  that  the  PN  code  chip  rate  remains  at  the  nominal  value 
of  10.23  Mbps.  This  assumption  is  important  for  the  interpretation  of  both 
experimental  and  simulation  results  for  receiver  performance  in  terms  of  actual 
GPS  performance  in  a  nuclear  environment.  An  unmodified  receiver  allows  for 
a  direct  study  of  GPS  X-sct  operation  in  a  nuclear  environment.  Of  course, 
if  receiver  modifications  wore  instituted  (for  example,  a  lower  code  chip 
rate),  the  numerical  simulations  could  be  easily  modified  accordingly. 

Step  1,  above,  involves  the  analysis  of  data  taken  from  three  air¬ 
craft  passes  in  the  shadow  of  two  barium  clouds,  Esther  and  Fern.  The 

* 

coherently  received  aircraft  data  was  then  numerically  back-propagated  to 
the  point  where  amplitude  scintillations  were  minimal  at  the  barium  cloud 
locations.  The  phase  of  the  signal  at  this  location  can  then  be  directly 
interpreted  as  electron-density  fluctuations  integrated  along  the  lino  of 
sight  through  the  barium  cloud.  This  phase  was  then  detrended  and  Fourier 
analyzed  to  obtain  measurements  of  integrated  electron-density  power  spectra, 
including  spectral  index,  phase  standard  deviation,  and  outer  scale. 

The  integrated  electron-density  power  spectral  density  (PSD)  is 
then  used  as  the  basis  of  a  statistical  model  of  barium  cloud  striations.  The 
model  includes  a  large  gaussian  shape  to  represent  the  mean  or  deterministic 
(unstriated)  portion  of  a  barium  cloud  plus  a  random  part  based,  in  part,  on 
the  measured  PSD  parameters  discussed  above.  The  inner  scale  is  included  as 

*  Dr.  C.  W.  Prettie  of  ESL,  Inc.,  provided  the  back-pi’opagation  results. 


an  additional  adjustable  parameter.  Inclusion  of  the  large  gaussian  profile 
to  represent  the  mean  barium  cloud,  while  providing  a  good  model  of  the  actual 
barium  propagation  environment,  causes  difficulty  in  interpretation  of  the 
propagation  results  and  related  receiver  simulation  results  because  of  the 
inhomogeneous  propagation  environment.  For  this  reason,  some  simulations  were 
made  without  the  deterministic  portion  of  the  barium  cloud  model.  In  Section  4 
propagation  results  are  presented  which  show  the  effect  of  both  the  deterministic 
and  random  portions  of  the  model. 

The  receiver  simulation  here  is  a  direct  digital  implementation  of 
the  X-set  receiver  and  includes  implementation  of  the  local  reference  generator/ 
correlator  and  signal  conditioner,  AFC/Costas  carrier  loop,  PN  code  loop,  AGC 
loop,  and  rate  multiplier/incremental  phase  modulators.  Also  included  are 
time-shared  PN  code  tracking  and  data  demodulation  and  decoding.  This  digital 
simulation  approach  provides  the  only  practical  method  of  accurately  analyzing 
the  operation  of  complex  nonlinear  receivers  in  the  presence  of  signal  propaga¬ 
tion  disturbances. 

No  effort  is  made  in  this  work  to  simulate  the  operation  of  the  GPS 
navigation  data  pi’ocessing  system.  Rather,  we  limit  consideration  to  the 
study  of  an  X-set  receiver  operating  in  a  disturbed  barium  environment— that 
is,  only  propagation  disturbances  on  a  single  channel  of  the  navigation 
processor  are  considered.  To  determine  the  effects  of  disturbed  propagation 
paths  on  the  GPS  navigation  processor,  our  X-set  receiver  model  would  have 
to  be  incorporated  into  a  larger  simulation  which  included  navigation  data 
processing. 

In  Section  4  of  this  report  several  detailed  examples  of  X-set 
receiver  operation  at  a  carrier  frequency  of  150  NBlz  are  pi'esented.  These 
examples  graphically  illustrate  the  phase,  frequency,  and  code  tracking 
performance  of  the  X-set  operating  in  a  striated  barium  environment. 
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Work  is  currently  in  progress  to  relate  propagation  results 
presented  in  this  report  to  analytical  results  for  mean  signal  arrival  time, 
time  delay  jitter  and  pulse  broadening  of  the  received  waveform.  These 
issues  will  be  addressed  in  a  separate  report. 
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SECTION  2 

BARIUM  CLOUD  STRIATION  MODEL 


The  multiple  phase  screen  (MPS)  propagation  simulation  provides  a 
numerical  solution  for  the  propagation  of  a  plane  wave  through  a  disturbed 
environment  characterized  by  index-of-refraction  fluctuations.  By  modeling 
a  barium  cloud  as  a  random  phase  screen,  it  is  possible  to  use  the  simula¬ 
tion  to  generate  signal  realizations  which  serve  as  direct  input  to  the  GPS 
X-set  receiver  simulation. 

This  section  describes  the  barium  cloud  striation  models  which 
were  developed  to  serve  as  input  to  the  MPS  propagation  simulation.  Results 
for  the  equivalent  phase  (integrated  electron-density)  power  spectral  den¬ 
sity  (PSD),  outer  scale,  and  spectral  index  are  presented.  Under  the 
assumption  of  a  10  km  propagation  distance  through  the  striated  barium 
cloud,  values  of  the  rms  electron-density  fluctuation  are  also  given. 

The  multiple  phase  screen  propagation  simulation  represents  the 
disturbed  region  by  a  number  of  phase  screens  located  between  the  satellite 
and  the  receiver.  Random  phase  fluctuations  in  each  screen  are  generated 
using  the  statistical  properties  of  the  electron-density  fluctuations  as 
determined  by  the  electron-density  power  spectral  density.  A  wave  (initially 
plane  as  it  enters  the  disturbed  region)  is  then  propagated  numerically 
from  one  screen  to  the  next  by  use  of  the  Fresnel-Kirchhoff  integral  equa¬ 
tion  until  a  solution  is  obtained  for  the  complex  electric  field  in  the 
receiver  plane.  This  technique  is  equivalent  to  a  solution  of  the  parabolic 
wave  equation  and  is  thus  able  to  account  for  multiple  scattering  (Reference 
2-1).  Since  the  phase  screens  are  random,  the  signal  propagated  to  the 
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receiver  is  random  and,  if  desired,  statistical  results  may  be  obtained  by 
averaging  a  number  of  different  simulations,  each  based  on  a  different 
sequence  of  random  numbers. 

The  random  phase  screens  are  chosen  in  a  statistical  manner  and 
are  characterized  by  a  phase  power  spectral  density  (PSD)  identical  to  the 
PSD  of  the  electron-density  fluctuations  integrated  along  the  direction  of 
propagation.  For  this  work,  the  phase  PSD  used  to  represent  a  striated 
barium  cloud  was  obtained  from  analvsis  of  observations  taken  during  the 
STRESS  barium  release  experiments  E;ther  and  Fern  (Reference  2-2). 


2.1 


BACK-PROPAGATED  PHASE 


This  section  describes  a  portion  of  the  propagation  experiments 
performed  during  the  STRESS  test  series  which  were  designed  to  obtain  in¬ 
formation  on  the  integrated  electron  density  fluctuations  produced  by 
barium  clouds.  Dr.  C.  W.  Prettie  of  ESL,  Incorporated  has  kindly  provided 
the  back-propagated  phase  data  discussed  below. 

In  order  to  observe  signal  scintillation  during  the  STRESS  experi¬ 
ment,  an  aircraft  was  flown  so  that  the  line-of-sight  from  the  aircraft  to 
the  geosynchronous  LES  8  satellite  intersected  the  striated  cloud.  The 
aircraft  velocity  was  perpendicular  to  the  projection  of  the  earth's  mag¬ 
netic  field  to  observe  the  small  scale  signal  fluctuations  caused  by  field- 
aligned  striations.  The  propagation  experiments  involved  either  uplink  or 
downlink  signals  at  different  frequencies.  In  the  uplink  experiment,  a 
340  MHz  CIV  signal  was  transmitted  from  the  aircraft,  up  through  the  barium 
cloud  to  LES  8,  then  was  relayed  along  an  undisturbed  path  to  a  rooftop 
receiving  station.  In  the  downlink  experiment,  a  250  MHz  CW  signal  was 
transmitted  from  the  rooftop  station,  up  to  LES  8,  then  down  thi'ough  the 
barium  cloud  to  the  aircraft.  For  both  the  uplink  and  the  downlink  experi¬ 
ments,  a  coherent  K-band  signal  was  used  as  a  phase  reference  so  that  the 
two  quadrature  components  of  the  lower  frequency  signals  could  be  measured 
and  recorded. 

During  the  three  aircraft  passes  of  interest  here,  the  received 
signal  exhibited  Rayleigh  fading  characteristics  with  many  deep  fades  and 
a  measured  scintillation  index  close  to  unity.  The  recorded  received 
signal  was  then  used  by  Dr.  Prettie  as  the  input  to  a  numerical  propagation 
simulation  where  the  signal  was  propagated  from  the  receiver  back  up  towards 
the  barium  cloud  (Reference  2-2).  This  "backwards"  propagation  in  free 
space  was  repeated  for  a  number  of  propagation  distances  until  a  minimum 
scintillation  index  was  obtained  for  the  back-propagated  electric  field. 


Under  ideal  circumstances,  where  the  barium-produced  phase  fluctuations  are 
perfectly  represented  as  a  thin  phase  screen,  then  the  fluctuations  of  the 
amplitude  of  the  electric  field  will  decrease  as  the  location  of  the  barium 
cloud  is  approached.  (Equivalently,  in  the  normal  forward  propagation  sense, 
a  thin  phase  screen  is  modeled  as  producing  instantaneous  phase  changes, 
but  no  amplitude  fluctuations.  The  amplitude  fluctuations  develop  as  the 
wave  propagates  away  from  the  phase  screen  location  (Reference  2-3)).  Now 
when  a  minimum  scintillation  index  is  obtained  for  a  particular  back- 
propagation  distance,  then  the  back-propagated  phase  at  this  distance  is 
assumed  to  be  directly  related  to  the  integrated  elects  on-density  fluctu¬ 
ations  caused  by  the  barium  cloud.  Reciprocity  permits  the  same  back- 
propagation  scheme  for  both  uplink  and  downlink  signals.  In  the  three 
cases  discussed  here  the  scintillation  index  obtained  after  back-propa¬ 
gation  is  approximately  0.5. 

The  usefulness  of  the  back-pi’opagation  technique  described  here 
is  dependent  on  a  number  of  factors.  First  of  all,  both  received  ampli¬ 
tude  and  phase  must  be  measured  sufficiently  accurately.  Secondly,  since 
the  numerical  propagation  simulation  hero  is  two-dimensional,  allowing 
for  variation  perpendicular  to  the  magnetic  field  but  permitting  no 
variation  along  the  magnetic  field,  the  technique  is  limited  in  the 
geometry  that  can  be  considered.  Fortunately,  the  experimental  geometry 
was  such  that  most  of  the  scattering  was  caused  by  irregularities  per¬ 
pendicular  to  the  magnetic  field  direction.  In  general,  however,  the 
effects  of  noise  and  finitr-s ized  correlation  lengths  along  the  magnetic 
field  have  not  been  investigated  here. 

Other  assumptions  implicit  in  this  analysis  are  that  the  barium 
cloud  is  relatively  stationary  during  the  aircraft  passes  and  that  very 
little  evolution  of  the  striations  occurs  over  the  pass  duration.  The 
first  assumption  is  reasonable  since  observations  of  cloud  motion  for  the 


two  clouds  considered  show  that  the  maximum  velocity  component  perpendicular 
to  the  projection  of  the  magnetic  field  line  on  the  earth's  surface  was 
approximately  18  m/sec  with  an  equal  velocity  component  in  the  orthogonal 
direction  along  the  earth's  surface.  The  second  assumption  also  seems 
reasonable  since  the  duration  of  the  three  aircraft  passes  are  70  seconds 
(Esther  pass  10  at  48  minutes  after  release),  130  seconds  (Esther  pass 
17  at  117  minutes  after  relase) ,  and  60  seconds  (Fern  pass  9  at  43  minutes 
after  release) . 


2.2 


PHASE  PSD  ANALYSIS 


This  section  describes  the  details  of  the  analysis  of  three 

measurements  of  back-propagated  phase  from  the  STRESS  series.  The  phase 

* 

measurements  were  provided  to  us  by  Dr.  Cliff  Prettie  of  ESL,  Incorporated 
(Reference  2-4).  As  stated  in  the  previous  section,  i";  is  assumed  that 
the  back-propagated  phase  is  directly  proportional  to  the  integrated 
electron-density  fluctuations  caused  by  the  barium  cloud  striations. 

Figure  2-1  shows  the  signal  phase  determined  from  the  ESL  back- 
propagation  technique  observed  during  an  aircraft  pass  in  the  shadow  of 
the  barium  cloud  Estlier  at  48  minutes  after  release. 

The  back-propagated  phase  data,  measured  as  a  function  of 
time  during  the  aircraft  pass,  was  converted  to  the  spatial  scale 
shown  on  Figure  2-1  by  obtaining  the  aircraft  velocity  perpendicular  to 
the  plane  containing  the  line-of-sight  and  the  magnetic  field  direction, 
both  measured  at  the  location  of  the  intersection  of  the  line-of-sight 
with  the  barium  cloud.  According  to  Reference  2-5,  if  the  striations  are 
a.xisymmetric  with  an  arbitrary  axial  ratio,  then  the  scale  size  perpendicular 
to  the  magnetic  field  direction  can  be  obtained  by  considering  scintilla¬ 
tions  only  along  this  particular  direction.  For  Esther  pass  10,  the  actual 
aircraft  velocity  was  235  m/sec  with  the  aircraft  travelling  at  a  constant 
altitude  along  a  direction  225°  clockwise  from  north.  The  aircraft  was 
located  at  29°N  latitude  and  85°K  longitude.  LES  8  was  in  a  synchronous 
orbit  at  21°iN  latitude  and  107°W  longitude  at  an  altitude  of  36000  km. 

At  the  barium  location  (180  km  altitude,  29°N  latitude,  86°K  longitude), 


Dr.  .f>rettie  also  provided  data  on  the  link  geometry  during  STRESS. 
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propagated  phase  for  Esther  pass  10. 
provided  by  ESL,  Incorporated.) 


the  magnetic  field  dip  angle  was  59°  and  the  angle  between  the  line-of-sight 
and  the  magnetic  field  direction  was  30°.  The  angle  between  the  aircraft 
direction  and  the  normal  to  the  plane  containing  the  line-of-sight  and  the 
magnetic  field  direction  was  67°.  By  design,  the  aircraft  velocity  Vv'as 
in  the  same  direction  as  the  horizontal  projection  of  the  vector  perpendi¬ 
cular  to  the  plane  containing  the  magnetic  field  direction  and  the  line-of- 
sight. 

The  pha^'e  data  is  easily  detrended  or  separated  into  a 
"deterministic"  or  non-striated  gaussian-like  part  and  a  "random"  or 
striated  part.  Figure  2-2  shows  the  deterministic  part  which  is  obtained 
from  the  data  of  Figure  2-1  by  consti'ucting  a  running  average  over  a 
3-kilometer  interval.  The  random  part  is  then  obtained  by  subtracting 
the  running  average  from  the  original  data.  This  procedure  of  separating 
out  the  detoiministic  part  of  the  phase  is  quite  reasonable  here  since 
aircraft  passes  at  earlier  times  after  release  quite  clearly  show  the 
smooth  gaussian-ltkc  unstriated  barium  cloud,  quite  similar  to  that  .>hown 
in  Figure  2-2. 

Once  the  random  part  the  phase  shown  in  Figure  2-i  wa.s  obtained, 

the  central  portion,  from  6  km  to  18  km  on  the  figure  was  then  iolected  as 

best  representing  the  homogeneous  or  unifoi’m  striation  region  of  the  cloud 

where  the  noticeably  different  appearance  of  the  edges  was  absent.  The 

2 

variance  of  the  phase,  was  obtained  from  this  portion  by  ..omputing  the 
necessary  averages  directly  from  the  data.  The  power  spectrum  shown  in 
Figure  2-4  was  also  obtained  from  this  unifoi-m  portion  of  the  data,  where 
the  "raw"  spectrum  was  smoothed  by  forming  a  running  average  over  a 
±  5%  bandwidth. 

Another  important  consideration  should  be  mentioned  regarding  the 
choice  of  the  uniform  striation  region  from  the  random  part  of  the  phase. 

In  addition  to  choosing  that  region  from  the  center  where  the  phase 
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fluctuations  appear  homogeneous,  it  is  also  convenient  to  choose  the  end 


r  points  carefully  so  that  no  discontinuity  exists.  Elimination  of  any 

discontinuity  at  this  stage  of  the  analysis  facilitates  a  cleaner  use 
-  of  a  digital  Fourier  transform  algorithm  to  compute  the  power  spectral 

!  density.  No  windowing  is  needed  and  any  difficulties  involving  rescaling 

1^  of  the  variance  because  of  windowing  effects  are  avoided.  As  a  check  on 

i  the  results,  the  PSD's  computed  from  the  carefully  chosen  uniform  regions 

•  were  compared  to  the  PSD's  of  the  original  data  before  removal  of  the 

deterministic  parts.  Agreement  was  noted  for  the  portion  of  the  PSD 
unaffected  by  the  detrending. 
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The  noise  level  noted  on  Figure  2-4  was  obtained  from  consulta¬ 
tions  with  Dr.  Prettie  of  ESL  (Refei'ence  2-4).  Also  shown  in  Figure  2-4 
is  the  best  fit  one-dimensional  power  spectral  density  (PSD)  of  the  form 
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where  is  the  variance,  obtained  as  explained  above.  T  is  the  gamma 
function,  is  the  outer  scale  and  K  is  the  spatial  wavenumber.  The 
form  of  the  PSD  is  chosen  so  that  the  normalization  is 
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The  fit  shown  was  obtained  by  determining  the  spectral  index  m  by  a  least 
square  fit  of  the  smoothed  PSD  between  K  =  lO"  and  10  radians/m.  In 
obtaining  the  spectral  index,  points  chosen  from  the  smoothed  PSD  were 
equally  spaced  on  the  wavenumber  axis  shown  so  that  the  results  were  not 
driven  by  the  large  number  of  data  i)oints  at  the  higher  wavenumber  end  of 
the  spectral  interval  over  which  the  least  squares  fit  was  performed. 
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Once  the  spectral  index,  m,  and  the  variance,  a^,  were  both  obtained,  the 
value  of  the  outer  scale  was  obtained  by  solving  Equation  2-1  at  K  =  0.1 
radians/m  with  S(K)  obtained  from  the  least  squares  fit. 

The  same  analysis  procedure  was  used  for  two  other  aircraft  passes 
during  releases  Esther  and  Fern.  Table  2-1  shows  the  results  for  the  spectral 
index  and  outer  scale  for  all  three  passes  analyzed.  As  stated  previously, 
the  phase  from  Esther  pass  10  was  detrended  by  subtractinj,  a  running  average 
obtained  over  a  3  km  interval.  For  Esther  pass  17  the  interval  chosen 
was  6  km;  for  Fern  pass  9,  the  appropriate  interval  was  2  km.  Figures  2-5 
and  2-6  show  the  PSP’s  of  the  random  parts  of  the  back-propagated  phase 
obtained  from  Esther  pass  17  and  Fern  pass  9.  The  link  geometry  for  Esther 
pass  17  and  Fern  pass  9  was  very  similar  to  the  geometry  of  Esther  pass  10 
and  similar  scaling  of  the  recorded  back-propagated  phase  was  used  to  obtain 
the  distance  scales. 


Table  2-1.  PSD  measurements  for  STRESS  events. 
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Note  that  none  of  the  three  PSD's  shown  in  Figures  2-4  to  2-6 
show  any  evidence  of  an  inner  scale  cutoff.  It  is  well  known  that  the 
inner  scale  plays  an  important  role  in  cutting  off  angular  scattering— 
particularly  for  phase  PSD's  with  a  spectral  index  of  3  or  less  (Ref¬ 
erences  2-6  and  2-7). 

It  should  be  noted  that  the  results  in  Table  2-1  include  the 
effects  of  a  careful  consideration  of  the  actual  experimental  geometry 
(aircraft  position  and  heading,  satellite  location,  magnetic  field  direction). 
Earlier  results,  based  on  less  detailed  information  on  the  experiment  and 
a  slightly  simpler  analysis  technique  (visual  instead  of  least  squares  curve 
fitting)  were  used  in  the  rest  of  this  report  to  provide  statistical  models 
of  the  important  features  of  Esther  pass  10  and  Fern  pass  9,  These  results 
differ  only  slightly  from  those  presented  in  Table  2-.',  For  Esther  pass  10 
the  original  value  of  was  360  m  with  a  spectral  index  of  3.  For  Fern 
pass  9,  the  original  value  obtained  for  was  270  meters  with  a  spectral 
index  of  3.5. 


Measurements  taken  during  the  STRESS  test  (Reference  2-8)  indicate 
that  the  thickness  of  the  striated  barium  clouds  range  from  5  to  10  kilometers 
perpendicular  to  the  magnetic  field  lines.  Of  course  elevation  angles  other 
than  90°  (vertical  propagation)  will  give  a  larger  propagation  distance 
through  the  striations.  For  simplicity,  assume  that  the  combined  effect 
of  cloud  thickness  and  elevation  angles  gives  a  propagation  distance  of 
10  km  through  the  barium  striations.  It  is  possible  to  deduce  the  mean- 
square  electron  density  fluctuations  from  the  relationship 

oj  .  2(r^X)2  ^  r  (2-3) 

where  r^  is  the  classical  electron  radius,  X  the  wavelength,  the 
outer  scale  size,  and  L  the  propagation  distance  through  the  striated 
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layer  (assumed  to  be  10  km).  Using  the  results  shown  in  Table  2-1  with  i 

Equation  2-3  the  values  obtained  for  the  rms  electron  density  fluctuations  I 

are  1.4x10^  cm  ^  for  Esther  pass  10,  9.8x10  cm"^  for  Esther  pass  17  and  j 

5-3  ^ 

7.0x10  cm  for  Fern  Pass  9,  These  values  are  within  the  ranges  ' 

1 

measured  by  radar  (Reference  2-9)  and  by  in-situ  rocket  probes  (Reference  ? 

2-10).  I 
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BARIUM  CLOUD  MODELING  VERIFICATION 


As  a  check  on  the  barium  cloud  striation  modeling,  early  results 
from  the  back-propagation  phase  analysis  of  Esther  pass  10  were  used  in  the 
multiple  phase  screen  simulation  to  generate  signals  to  compare  to  the  actual 
received  signals  measured  experimentally.  These  early  results  differed 
from  those  in  Table  2-1  in  that  the  slopes  were  obtained  by  a  visual  fit 
to  the  PSD's.  For  Esther  pass  10  the  early  value  of  the  outer  scale  was 
360  III  (as  compared  to  390  m  obtained  by  the  least  squares  process) .  The 
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spectral  index  of  K  did  not  change  with  the  more  accurate  analysis. 


Figure  2-7  shows  the  phase  screen  generated  to  model  Esther  pass  10 
and  is  to  bo  compared  to  the  back-propagated  phase  data  shown  in  Figure  2-1. 
The  phase  screen  shown  in  Figure  2-7  is  composed  partly  of  a  gaussian  shape 
whose  size  is  chosen  to  match  the  deterministic  barium  cloud.  To  this 
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deterministic  gaussian  is  added  a  random  phase  characterized  by  a  K  PSD 
with  an  outer  scale  of  360  m.  The  random  phase  is  windowed  or  cutoff  on 
the  edges  so  that  no  phaso  pertur’bations  occur  outside  the  deterministic 
gaussian  cloud.  A  pseudo-random  sequence  of  coefficients  actually  determines 
the  details  of  the  random  phase  screen  so  that  the  phase  screen  realization 
shown  in  Figure  2-7  is  one  of  many  possible.  To  obtain  statistically  sig¬ 
nificant  results,  other  pseudo-random  sequences  could  be  used  and  the  results 
for  the  received  signal  could  be  averaged. 

To  model  Esther  pass  10,  a  plane  wave  is  propagated  through  the 
phase  screen,  down  to  a  ground  based  receiver.  16384  points  were  used  to 
represent  the  phase  shown  in  Figure  2-7  over  a  distance  of  30  km,  of  which 
only  20  km  is  shown  in  the  figure.  Thus,  the  distance  between  points  on 
the  grid  is  1.831  m.  In  the  absence  of  an  explicit  inner  scale  in  the  PSD, 
the  minimum  separation  between  grid  points  controls  the  effective  inner  scale 
in  the  MPS  code.  Here  the  effective  inner  scale  is  Ax/tt  or  0,58  m,  accord¬ 
ing  to  the  Nyquist  sampling  theorem. 
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Figures  2-8  and  2-9  show  a  comparison  of  the  numerically  generated 
signal  amplitude  (plotted  in  decibels  relative  to  the  mean  amplitude)  to  the 
actual  experimentally  received  data.  Figure  2-8  shows  23  seconds  (4.6  km 
at  a  velocity  of  220  m/sec)  of  received  data  taken  during  the  STRESS 
test.  Figure  2-9  shows  a  3  km  portion  of  the  signal  realization  generated 
numerically.  The  15  km  point  on  the  abscissa  of  Figure  2-9  corresponds 
to  the  point  directly  below  v.he  center  of  the  cloud.  A  visual  comparison 
of  the  two  figures  indicates  x'easonably  good  agreement  of  the  numerical 
results  with  the  observed  experimental  results.  The  major  difference 
appears  to  be  that  the  numerical  results  have  a  higher  percentage  of  fades 
below  20  dB.  The  distance  between  fades  appears  similar.  No  statistical 
comparison  was  possible  here  since  the  experimental  data  was  available 
only  in  the  graphical  form  shown. 

Figures  2-10  and  2-11  compare  the  numerical  and  measured  signal 
phases  corresponding  to  the  previous  signal  amplitudes.  Note  a  slight 
discrepancy  in  the  ordinate  axes  where  t.ie  two  scales  differ  in  length  by 
2Tr  radians.  A  visual  comparison  again  indicates  relatively  good  agreement 
between  overall  appearances  of  the  numerical  and  experimental  results. 

Figures  2-12  and  2-13  compare  the  experimental  back-propagated 
★ 

phase  for  Fern  pass  9  with  the  numei'ical  phase  screen  model  of  this  pass. 
Figure  2-12  shows  the  back-propagated  phasi  data  provided  by  Dr.  Prettie 
of  ESL.  By  a  visual  fitting  process  an  outer  scale  of  270  meters  and  a 
spectral  index  of  3.5  were  originally  obtained  for  the  phase  PSD.  (More 
detailed  analysis  indicates  that  a  better  fit  is  obtained  for  an  outer 
scale  of  350  m  and  a  spectral  index  of  3.2  as  shown  in  Table  ^-1.) 

*  A  velocity  of  220  m/sec  was  originally  used  to  convert  the  tine  domain 
data  of  Fern  pass  9  to  the  spatial  domain.  More  careful  attention  to 
the  actual  geometry  later  resulted  in  a  value  of  200  m/sec  which  was 
used  to  obtain  the  more  accurate  results  shown  in  Table  2-1. 
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Experimentany  measured  received  signal  amplitude  during  Esther 
(Figure  provided  fay  ESL,  Incorporated.) 


Experimentally  measured  signal  phase  corresponding  to  signal  amplitude  shown 
in  Figure  2-8.  (Figure  provided  by  ESL,  Incorporated.) 


Again  using  16384  points  to  represent  a  grid  30  km  long,  the 
original  PSD  parameters  were  used  to  generate  the  phase  screen  shown  in 
Figure  2-13  as  a  match  to  the  experimental  data.  A  comparison  of  Figures 
2-12  and  2-13  demonstrates  the  reasonableness  of  the  original  fit 
to  the  PSD. 


Figures  2-14  and  2-15  compare  the  experimentally  measured  received 
signal  amplitude  for  Fern  pass  9  witli  the  amplitude  obtained  after  numer¬ 
ical  propagation  through  the  phase  s'reon  shov^^n  in  Figure  2-13.  A  visual 
comparison  of  Figures  2-14  and  2~15  agoin  show  the  reasonableness  of  the 
phase  PSD  analysis  and  of  the  numerical  propagation  procedure. 

Figures  2-16  and  2-17  comarro  the  corresponding  experimentally 
measured  received  phase  with  the  numerically  determined  received  phase. 
Comparison  of  these  two  figures  also  indicates  reasonably  good  agreement 
between  data  and  simulated  signal  characteristics. 

Thus,  this  section  demonstrates  the  general  adequacy  of  the 
results  of  the  back-propagated  phase  analysis.  When  the  back-propagated 
phase  PSD  is  used  as  input  to  the  MPS  propagation  simulation,  the  results 
for  the  received  signal  amplitude  and  phase  are  in  good  agreement  with  the 
experimental  results.  With  the  e.xception  of  the  inner  scale  values  and 
the  actual  propagation  geometry.  Table  2-1  lists  all  the  necessary  para¬ 
meters  to  utilize  the  MPS  propagation  simulation  to  generate  signal 
structures  characteristic  of  propagation  through  a  barium  cloud  for  the 
PLACES  experiment. 

In  Section  4,  the  original  results  for  the  back-propagated  phase 
PSD  analysis  of  Esther  pass  10  and  Fern  pass  9  arc  used  to  generate  signal 
structures  to  represent  propagation  through  the  barium  clouds  at  carrier 
frequencies  from  150  to  250  NDIz. 
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TIME  (1  SECOND  PER  TICK) 

Experimentally  measured  received  signal  amplitude  during  Fern  pass 
(Figure  provided  by  ESL,  Incorporated.) 


mentally  measured  signal  phase  corresponding  to  signal  amplitude  shown 
lure  2-14.  (Figure  provided  by  ESL,  Incorporated.) 


Figure  2-17.  Numerically  generated  signal  phase  corresponding  to  signal  amplitude  shown  in 
Figure  2-15. 
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SECTION  3 

PN  SPREAD  SPECTRUM  CODE  CORRELATOR 


In  this  section  a  technique  is  derived  to  accurately  simulate  the 
operation  of  the  pseudonoise  (PN)  code  correlator  in  the  X-set  spread  spectrum 
receiver.  With  the  addition  of  the  algorithm  presented  here,  the  X-set 
simulation  described  in  Reference  3-1  is  capable  of  fully  simulating  X-set 
receiver  performance  in  all  types  of  conditions,  including  severe,  frequency- 
selective  environments.  In  the  next  section,  frequency-selective  multiple 
phase  screen  (MPS)  calculations  of  signal  propagation  through  barium  clouds 
are  used  as  direct  input  to  the  X-set  receiver  simulation  and  results  pre¬ 
sented  for  receiver  performance. 

Frequency-selective  effects  can  degrade  the  performance  of  the 
PN  code  tracking  loop  in  a  spread  spectrum  receiver.  Signal  transmission 
through  a  striated  region  can  cause  time  delay,  time  delay  jitter,  and 
distortion  of  the  received  waveform,  all  of  which  affect  the  code  tracking 
performance.  Before  describing  the  technique  to  model  the  X-set  code  cor¬ 
relator,  it  is  useful  to  discuss  receiver  code  tracking  in  a  frequency- 
selective  environment. 
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3.1 


PN  CODE  TRACKING 


The  X-set  receiver  employs  a  noncoherent  delay-lock  loop  for 
pseudonoiss  code  tracking.  In  the  code  loop  the  incoming  signal  plus 
noise  is  alternately  correlated  with  early  and  late  codes  from  the  receiver 
code  generator.  This  is  known  as  tau-dithering,  with  the  amount  of  time 
dither  set  at  one-half  the  chip  duration  in  this  case.  The  dithering  loop 
is  basically  a  delay-lock  loop  in  which  the  early  and  late  local  references 
share  the  correlator  and  detector. 

The  receiver  estimates  of  time-delay  error  are  obtained  by  sub¬ 
tracting  measurements  of  late  power  (obtained  by  correlation  of  the  late 
code  with  the  received  code)  from  measurements  of  early  power  (obtained 
by  correlation  of  the  early  code  with  the  received  code). 

For  infinite  bandwidth,  and  no  scintillation,  dispersion  or  noise, 
the  envelope  of  the  PN  code  correlation  function  formed  by  correlating  the 
received  code  with  the  internal  reference  code  is  a  perfect  triangle  as  shown 
in  Figure  5-1.  For  perfect  code  tracking,  the  ^rly  power  is  obtained  from 
Figure  3-1  at  a  time  advance  of  one-half  chip  (v  =  -0.5)  and  the  late  power 
Is  obtained  from  the  figure  at  a  time  delay  of  one-half  chip  (x  =  0.5). 

Since  these  values  are  identical,  subtracting  the  late  power  from  the  early 
power  yields  zero  time-delay  eri'or,  as  it  should  for  perfect  code  tracking. 


-1  0  1 
TIME  DELAY  (CHIPS) 

Figure  3-1.  Envelope  of  the  PN  code  autocorrelation  function,  infinite  band¬ 
width,  no  scintillation  or  dispersion. 
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Figure  3-2  shows  a  three-dimensional  plot  of  the  PN  code  autocorre¬ 
lation  function  envelope  for  nonselective  fading  conditions.  That  is,  all  fre¬ 
quency  components  across  the  bandwidth  of  the  PN  spread  spectrum  signal  fade 
identically  in  this  example.  The  PN  code  chip  rate  is  :r.uch  larger  than  the 
fading  rate  here  so  that  propagation  conditions  remain  unchanged  during  the 
many  chip  durations  required  to  calculate  one  of  the  PN  code  autocorrelation 
functions  shown.  As  time  progresses,  fading  occurs  and  subsequent  measure¬ 
ments  of  the  autocorrelation  function  are  shown  plotted  behind  previous 
measurements  in  the  figure.  For  nonselective  fading,  the  shape  of  the 
autocorrelation  function  remains  triangular  although  the  peak  and  sharp 
edges  are  rounded  because  of  the  finite  signal  bandwidth.  Fading  effects 
are  shown  in  Figure  3-2,  by  noting  the  progression  of  the  envelope  peaks 
witn  time. 

Frequency-selective  fading  conditions  caused  by  severe  electron- 
density  fluctuation^'  can  cause  notable  time  delay  jitter  for  wide  bandwidth 
systems  in  two  ways.  First,  striations  can  cause  frequency-dependent  angular 
scattering  where  stronger  scattering  occurs  for  the  lower  frequency  components. 
Secondly,  during  propagation  through  an  ionited  medium,  lower  frequency  signal 
components  will  experience  a  greater  time  delay  than  higher  frequency  compo¬ 
nents  . 

Both  frequency-selective  propagation  mechanisms  contribute  tc  the 
results  in  Figure  3-3  which  shows  the  PN  code  autocorrelation  function 
behavior  over  a  time  interval  of  several  seconds.  As  can  be  seen  in  this 
example,  the  PN  code  correlation  function  is  severely  distorted  with 
signal  energy  delayed  up  to  four  code  chips. 

Figure  3-3  also  shows  a  rather  unusual,  although  interesting, 
example  of  this  severe,  frequency-selective  environment  on  the  operation 
of  the  PN  code-tracking  function.  The  arrows  in  Figure  3-3  show  the 
estimated  code  time  delay  as  it  is  tracked  in  the  receiver  delay- lock  loop. 
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Figure  3-3.  PN  code  autocorrelation  function  in  a  frequency-selective 

environment.  The  arrows  show  the  locus  of  the  peak  tracked 
by  the  receiver  code  loop. 
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Near  the  front  of  the  figure,  the  time  delay  estimate  is  about  one-half  chip 
with  the  receiver  code  loop  correctly  tracking  the  peak  of  the  code  correlator 
output.  However,  in  this  example,  the  details  of  the  code  correlation 
function  time  development  cause  the  code  loop  to  follow  a  misleading  path, 
indicated  by  the  arrows  in  the  three-dimensional  plot,  so  that  the  time 
delay  error  increases  to  about  two  chips  during  tlie  brief  time  interval 
shown.  At  this  point,  the  code  loop  is  no  longer  tracking  the  main  peak  of 
the  code  correlation  function,  but  rather  it  is  following  an  extended  side- 
lobe.  (Later  in  this  particular  simulation  the  delay  estimate  drops  back  to 
around  one-half  chip.) 


As  can  be  seen  from  the  discussion  above,  accurate  calculation  of 
the  P.N  code  correlator  output  is  necessary  to  adequately  simulate  the 
operation  of  the  X-set  receiver.  In  the  following,  a  method  is  given  to 
realistically  model  tlie  X-set  P.N  code  correlation  operation  in  the  presence 
of  frequency-select i ve  scintillation  and  dispersion. 


3.2 


THEORETICAL  DEVELOPMENT  -  RECEIVED  WAVEFORM 


Let  the  real  waveform  transmitted  by  the  satellite  be  expressed 


as 


x(t) 


(  t) 

=  Re  I  m(t)e  °  j 


I 

1  ,,,  jV  1 
=  2 


+  J  m  (t)e 


(3-1) 


where  oj^  Is  the  carrier  angular  frequency  and  m(t)  is  the  modulation 
waveform.  For  a  pseudonoise  (PN)  code,  m(t)  consists  of  a  long  sequence 
of  pseudorandom  polarity  changes  (binary  PSK-modulated  signal).  The 
frequency-domain  representation  of  the  transmitted  waveform  is  given  by 


X(co) 
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where  M(w)  is  the  Fourier  translarm  of  the  modulation  waveform: 

V 

M(uj)  =  J"  m(t)e  dt  (3-3) 

— >Jt> 

.Now  after  propagation  through  an  loniccd  medium  for  a  distance 
z  ,  the  signal  spectrum  is  given  by 


Vij)  =  \(u)e’ 


(3-4) 
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where  the  term  exp[j6((ji))]  is  the  transfer  function  of  the  ionized  medium. 
For  a  "smooth"  plasma,  with  no  striations  and  negligible  electron  collisions, 


the  transfer  function  is 


2  1/2 


where  is  the  propagation  distance  through  the  ionized  medium,  c  is 

the  speed  of  light,  and  is  the  plasma  frequency.  In  terms  of  the  mean 

electron  density  N 


9  9 

0)“  =  4'rTc‘r  N 
])  c  e 


(3-6) 


wl\ere  r  is  the  classical  electi’on  radius  (2.82x10  m) .  Thus  after 
e 

propagation  through  a  smooth  plasma,  the  spectrum  of  the  received  waveform 
is  obtained  from  I-quations  3-2  and  3-4  as 


'(to)  =  ^  M(a)-a)^)e-^*^^*^^  +  y  M  (-w-a)^)e 
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Now  in  the  first  term  of  liquation  3-7  expand  0((Ji))  in  a  Taylor 


series  about  o)  : 
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whore  the  subscripted  plus  sign  signifies  the  expansion  about  a  positive 
carrier  frequency  and  the  primes  denote  differentiation  with  respect  to  to. 
From  equation  3-5 
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In  the  second  term  of  liquation  3-7,  expand  0(o))  in  a  Taylor  series  about 
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Por  notational  simplicity,  let  0(a)^)  =  0^,  0  (u)^)  =  Qq  ^o  forth. 

Retaining  tl)c  first  four  terms  in  the  Taylor  series,  the  received  wave¬ 
form  for  smootli  plasma  effects  is 
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do) 
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In  the  first  integral  use  the  substitution  v  =  w-w^.  In  the  second  integral 


let  V  =  -co-u^.  Then  Equation  3-14  becomes 
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where  0^  can  be  seen  to  be  the  total  mean  phase  shift  at  the  carrier 


frequency  and  t^  is  the  mean  group  time  delay  (i.e.,  the  group  delay 


at  the  carrier  frequency): 
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l-'or  cases  of  interest  .>nere  ic“  ‘^p  (otherwise  the  wave  cannot  propagate 


witliout  severe  attenuation),  0  and  t  can  be  exiiressed  as: 
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where  the  first  terms  above  are  the  frce-space  contribution,  and  the 
>econd  terms  are  jiroport ional  to  the  total  electron  concent.  Note  that 
the  second  term  in  liiuation  3-15,  arising  from  the  negative  frequenev 
spectral  components,  is  the  comple.x  conjugate  of  the  first  term,  'ilius 
the  received  waveform  v(t)  i'^  real,  as  it  should  be. 
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Now  consider  the  possible  frequency-selective  effects  of  a  stri¬ 
ated  plasma.  Let  H(v),  where  v  =  a)-a)^,  denote  the  complex  signal  received 
after  propagation  of  a  plane  wave  at  a  positive  frequency  oj  through 
the  striated  plasma.  As  illustrated  in  Reference  3-2,  H(v)  is  the  received 
signal  computed  from  a  multiple  phase  screen  (MPS)  calculation  of  propaga¬ 
tion  through  an  ionized  region.  It  is  clear  from  Equation  3-15  that  H(v) 

•k 

modifies  the  positive  frequency  spectral  components  and  that  H  (v)  modifies 
the  negative  frequency  spectral  components.  Thus  the  received  waveform 
including  both  smooth  and  striated  plasma  disturbances  is  given  by 
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This  waveform  at  the  receiver  input  can  be  written  in  terms  of 
a  complex  modulation  envelope, 
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where  the  complex  modulation  envelope  is 


(3-20) 
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Note  that  Lquation  3-20  has  the  same  general  form  as  Equation  3-1  —  a  wave 
modulated  on  a  carrier  at  angular  frequency  witli  an  additional  time  delay 

and  phase  shift  as  given  by  Equations  3-17  and  3-18. 
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3.3 


X-SET  PN  CODE  CORRELATION  OPERATION 


The  X-set  receiver  simulation  documented  in  Reference  3-1  is  a 
direct  sampled-data  digital  implementation  of  the  GPS  X-set  receiver 
described  in  Reference  ,3-3.  The  simulation  includes  implementation  of  the 
local  reference  generator/correlator  and  signal  conditioner,  AFC/Costas 
carrier  loop,  PN  code  loop,  AGC  loop,  and  rate  multiplier/incremental  phase 
modulators.  Figures  3-4  and  3-5  show  functional  diagrams  of  the  implementa¬ 
tion  of  the  AFC/Costas  loop  and  the  PN  code  tracking  loop,  respectively. 

The  details  of  the  various  frequency  conversions  preceding  the 
inphase  and  quadrature  channel  product  detectors  are  treated  in  Reference  3-1 
but  are  unimportant  to  the  analysis  hero.  For  the  following  development,  the 
simple  block  diagram  of  Figure  3-6  will  suffice. 


In  the  X-set  spread  spectrum  receiver,  the  received  waveform 
v(t)  is  multiplied  by  a  time-shifted  replica  of  the  PN  modulation  waveform 
m(t),  which  has  been  modulated  onto  a  local  reference  carrier.  The 
locally  generated  real  PN  wavefoi-m  is  represented  by 

j(a)  t+0  (t))  *  -j(w  t-iO  (t)) 

u^(t)  =  ni(t-tj)  e  +  111  (t-t^)  e  (3-22) 

where  t^  is  the  receive*’  estimate  of  the  propagation  time  delay,  0j^(,t) 
is  the  receiver  estimate  of  signal  phase  shift,  and  is  the  local 

oscillator  freciuency  co^) . 
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'K  f  i^' 


where  the  sum-frequency  components  have  been  removed  by  filtering,  leaving 
only  the  difference-frequency  components.  The  difference  frequency  is 
given  by 

(i).  =  w  -u),  (3-24) 

2  0  1 

and  (p  is  the  error  in  the  receiver  estimate  of  carrier  phase: 


4.  =  <})(t)  =  0^(t)  -  e^(t) 


(3-25) 


The  waveform  w(t)  is  now  mixed  with  two  quadrature  sinusoids, 
cos{i^2^  -sinw2t,  to  form  two  baseband  channels  in  the  receiver.  After 

double-frequency  components  are  filtered  out,  these  baseband  voltages  are 
given  by: 


if  *  * 

i(t)  =  2  ® 


-tj)  e  J 


if  *  J*)*  * 

q(t)  =  •^j^E(t-t^)  m  (t-tj)  e  -  E  (t-t^)  m(t-tp  e  J 


(3-26) 


(3-27) 


These  two  quadrature  channel  signals  can  be  rewritten  as 


i(t)  =  Re  I  E(t-t^)  m  (t-t^  c 


(3-28) 


q(t)  =  Im  j  E(t-t^)  m  (t-t^)  e 


(3-29) 


The  two  baseband  voltages  are  integrated  over  a  time  interval 
At  to  foi’m  I  and  Q  samples  for  subsequent  digital  processing  in  the 
receiver.  Since  integration  is  a  linear  process  it  is  clear  that  the  I 
and  Q  samples  are  given  by  the  real  and  imaginary  parts  of  the  integral: 
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^«ig»rin  I'wrntsa’ifff^j 
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'  *  '  id) 

E(t  +t)  m  (t  )  dt' 


(3-30) 


4^. 

2 


where  t  =  t-t^  references  the  time  to  the  receiver,  and  t  is  the  error 


in  the  receiver  estimate  cf  propagation  time  delay  through  the  medium, 


T  =  t,-t  .  As  stated  earlier,  t  is  the  satellite  or  transmitter  time, 
1  0  ’  ’ 


tj^  is  the  receiver  estimate  of  propagation  time  delay,  and  t^  is  the 
actual  propagation  time  delay. 


Using  Equation  3-21  in  Equation  3-30,  the  equation  for  p 


becomes 


4i 

?  r 


P(T) 


•/■  */ 


45 

9 


M(v)  l'(v)  e 


jv(t  +T) 


dv 


m  (t  ie"^  dt 


(3-31) 


where  F(v)  is  the  total  plasma  frequency-selective  transfer  function, 
including  both  smooth  and  striated  plasma  components: 


F(v)  =  ll(v)  e 


1  II  9  1 

j  (4  0  v“  +  f  e  V  ) 

•'^2  o  0  0  ^ 


(3-32) 


if  the  propagation  disturbances  do  not  vary  significantly  with  time  over 
the  receiver  sampling  period,  that  is,  if  F(v)  and  (p  are  essentially 
constant  over  the  time  interval  At,  then  Equation  3-32  can  be  rewritten 
as 
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MC")  F(v) 


r  ^ 

I  *  I 

J  m  (t  ) 


dt  dv 


(3-33) 


In  the  X-set,  the  integration  period  At  is  4  ms  which  en¬ 
compasses  40,920  chips  of  the  PN  code  sequence.  In  this  case,  the  limits 
on  the  t'  integral  can  be  set  to  infinity  with  little  error.  Thus 


M*(v)  =  lim  r  m*(t’)e^^^  dt' 

At-x»  J 


(3-34) 


and  Equation  3-33  reduces  to 


*  ^ 

=  V  /  |M(v)|“  F(v)  e-'^Vj 


(3-35) 


|M(v)|  ,  the  squared  magnitude  of  the  Fourier  transform  of  the  modulation 

waveform,  is  the  power  spectrum  of  the  PN  sequence,  given  by  the  Fourier 

transform  of  the  PN  code  autocorrelation  function  (1  -  |t|/T  ),  where  T 

‘  •  c  c 

is  the  PN  code  chip  period  (Reference  3-4,  p.  341),  Thus 


/  sinvT  /2  \  2 


•’'cl -rr} 


(3-36) 


As  a  check  for  nonselective  cases  (no  scintillation  and  no 
dispersion),  F(v)  is  unity,  and  the  use  of  Equation  3-36  in  Equation 
3-35  yields. 


P(T)  .  -  M] 


(3-37) 


as  it  should. 
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In  general,  F(v)  '  not  unity  because  of  fading,  dispersion  and 
frequency-selective  scintil’  on  disturbances.  Equation  3-35  is  a  genera] 
expression  that  provides  the  PN  code  correlation  output  amplitude  and  phase 
(I  and  Q  baseband  channels)  for  any  type  of  propagation  condition. 
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COMPUTATIONAL  CONSIDERATIONS 


Equation  3-35  can  be  evaluated  by  discrete  Fourier  transform 
techniques.  The  X-set  receiver  requires  the  value  of  the  code  correlator 
output  at  only  one  (or  two)  specified  values  of  time  delay  (or  advance)  at 
a  given  simulation  time  (Reference  3-1).  Thus  the  use  of  a  simple  discrete 
Fourier  transform  (DFT)  algorithm  is  more  efficient  than  the  use  of  fast 
Fourier  transform  (FFT)  techniques.  Since  DFT's  do  not  assume  a  fixed 
sampling  interval  (as  do  FFT's),  this  procedure  is  more  general  than  FFT 
techniques  and  requires  less  CPU  time  for  this  particular  application. 


At  each  simulation  time,  an  array  of  Fourier  coefficients  is 
formed  at  a  sot  of  N  discrete  frequencies  across  the  PN  signal  bandwidth. 
Each  Fourier  coefficient  consists  of  three  multiplicative  factors:  the  PN 
signal  power  spectrum  |M(v)|  computed  from  Equation  3-36,  frequency- 
selective  propagation  factors  ll(v)  computed  using  the  MPS  technique,  anu 
dispersive  effects  duo  to  the  mean  total  electron  content  (TEC)  computed 
from  Equations  3-9  through  3-12.  The  phase  shift  and  time  delay  due  to  the 
mean  TEC  are  given  in  Equations  3-17  and  3-18.  The  corresponding  higher- 
order  dispersive  terms  arc-  found  from  Equations  3-11  and  3-12  to  be 
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where  it  has  again  been  assumed  that 
of  interest  here. 


9  9 

uj  “  »  0)  which  is  valid 
0  p  ’ 


in  all  cases 


Once  the  array  of  frequency-selective  Fourier  coefficients  C(n) 
has  been  computed  at  a  given  simulation  rime,  the  following  DFT  algorithm 
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is  used  to  compute  the  PN  code  correlator  output  amplitude  and  phase  at  any 
specified  value  of  time  delay: 

N/2-1  .  . . 

p(T)  =  Y.  +  C(N/2)cos(NTrAfT)  .. 

n=  -N/2+1 

where  Af  is  the  frequency  spacing  between  each  of  the  N  frequencies 
(N  even).  The  corresponding  unambiguous  range  of  time  delay  T  (DFT  compu¬ 
tational  window)  is  equal  to  1/Af. 

Frequency-selective  propagation  disturbances  introduce  additional 

mean  time  delay,  in  excess  of  that  given  by  Equation  3-18,  and  cause  p(t) 

to  be  generally  asymmetrical  and  shifted  toward  positive  t  (increased 

delay) .  Therefore  the  DFT  computational  window  is  not  centered  about  x  =  0 

but  rather  is  shifted  to  larger  values  of  delay  and  centered  around  a 

positive  value  of  x.  For  example,  with  32  frequencies  in  a  bandwidth  of 

2/T^,  where  T^  is  the  PN  chip  period,  the  DFT  window  is  set 

-2T  ^x^  14  T  . 
c  c 
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SECTION  4 

PROPAGATION  AND  X-SET  PERFORMANCE 


In  this  section,  examples  of  signal  scintillation  are  presented 
for  propagation  through  barium  clouds  at  frequencies  of  250,  200,  and  150  MHz. 
These  propagation  calculations  are  based  on  the  barium  cloud  models  discussed 
in  Section  2  and  represent  the  propagation  environment  for  two  different 
barium  clouds  (F.sther  and  I'em)  observed  at  48  and  43  minutes  after  barium 
release,  respectively. 

Detailed  examples  of  propagation  conditions  are  given  which  show 
the  development  of  the  received  signal  with  propagation  distance,  the  effect 
of  the  deterministic  barium  cloud,  and  the  effects  of  inner  scale  size. 
Examples  are  presented  showing  the  distortion  of  the  PN  code  correlator 
output  for  carrier  frequencies  of  150  and  200  MHz. 

This  section  also  includes  several  detailed  examples  of  the  per¬ 
formance  of  the  X-set  receiver  operating  at  a  carrier  frequency  of  150  MHz 
for  two  values  of  the  mean  carrier  power-to-noi se  density  ratio.  The 
examples  graphically  illustrate  the  phase,  frequency,  and  code  tracking 
perfomance  of  the  X-set  in  a  striated  barium  environment. 
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4.1 


NUMERICAL  BARIUM  CLOUD  .MODEL 


As  discussed  in  Section  2,  a  statistical  phase  screen  model  for 
barium  clouds  Esther  (48  minutes  after  release)  and  Fern  (43  minutes  after 
release)  was  developed  to  use  as  input  to  the  px'opagation  simulation.  The 
model  consists  of  a  deterministic  mean  barium  cloud  represented  by  a  large 
gaussian  shape  to  which  is  added  random  phase  irregularities  generated  using 
previously  developed  algorithms  (Reference  4-1).  The  random  phase  is  cutoff 
by  a  smooth  window  so  that  no  phase  irregularities  appear  outside  the  large 
gaussian  phase  structure.  This  phase  screen  representation  for  a  barium 
cloud  is  described  by  the  equation 


(p  e 
0 


-(x'.x  )^/r  2 
in  o 


-(x-x  +«.  )2/r.2 

+  (^(x)e  ^  U(.Xjj^-X.^-x) 


•t  <j>(x)e 


-(x-x  -£  )2/r  2 
m  c-^  c 


U(x-x  -i  ) 
^  m  c 


+  (|)(x)  [U(x-.x^^+£^) 


U(x 


•X  -I  )] 


(4-1) 


1  ,  X  ^  0 

U(x)  =  (4-2) 

(O  ,  X  <  0 

where  the  first  term  represents  the  deterministic  barium  cloud  with  a 
maximum  phase  of  (f  and  a  half-width  of  r  .  The  barium  cloud  is  centered 
at  x^^^  which  is  always  chosen  as  the  middle  of  the  MRS  grid.  ({i(x)  is  a 
random  phase  structure  which  is  generated  using  the  statistical  properties 
of  the  phase  fluctuations  as  determined  in  the  analysis  described  in 
Section  2.  Since  the  propagation  algorithm  generates  (l!(x)  as  a  homogeneous 
function  over  the  entire  MPS  grid,  simply  adding  ({)(x)  to  the  first  term 
in  Equation  4-1  would  yield  a  structural  barium  cloud  with  striations  off 
to  the  side.  Thus  the  last  three  terms  in  Equation  4-1  are  included  to  add 
a  windowed  version  of  (j)(x)  to  the  mean  gaussian  phase  representing  the 
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deterministic  barium  cloud.  As  is  evident  from  Equation  4-1  the  window  is 

unity  from  x  -  I  to  x  +  I  ,  and  is  cutoff  by  gaussian  functions  of  half- 
m  c  me  ^ 

width  r  beginning  at  x=x  -  £  and  x=x  +£  .  Figures  2-7  and  2-13 
c  ^'’  me  mc^ 

show  examples  of  the  phase  screen  models  of  Esther  at  48  minutes  after 
barium  release  and  Fern  at  43  minutes  after  barium  release.  Table  4-1  lists 
the  parameters  necessary  to  generate  phase  screen  models  of  barium  clouds 
with  the  MPS  propagation  simulation. 


Table  4-1.  Barium  cloud  parameters  used  in  predictions  for  the 
PLACES  experiment. 


Barium  Cloud 

- 1; — 

Lo 

m 

* 

d) 

^0 

'"o 

£ 

c 

r 

c 

Esther  (48  min) 

9.7  rad 

360  m 

3 

82  rad 

5  km 

5  km 

1  km 

Fern  (43  min) 

4.9  rad 

270  m 

3.5 

125  rad 

3  km 

4  km 

0.3  km 

*  0^  and  as  listed  correspond  to  a  frequency  of  340  MHz. 


As  discussed  above,  the  phase  fluctuations  are  represented  by  a 
power- law  PSD  of  the  form 


m  2  2 

o~  r(^)L^exp(-K"£.  ) 

/it  r(!^)(i+rL^“ 


(4-3) 


where  F  is  the  gamma  function,  K  is  the  spatial  wavenumber,  m  is  the 
spectral  index  and  is  the  outer  scale.  The  inner  scale,  ,  is  here 

used  explicitly  to  cutoff  the  spectrum.  In  some  cases,  the  value  is 

set  to  zero.  Under  this  circumstance,  the  value  of  the  inner  scale  is 
effectively  set  by  the  minimum  distance  between  points  on  the  MPS  grid.  As 
discussed  in  Section  2.3,  this  minimum  distance  of  1.831  m  corresponds  to 
an  effective  inner  scale  of  0.58  m. 
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We  have  found  that  utilization  of  the  complete  barium  cloud  model, 
including  the  mean  barium  cloud,  although  more  appropriate  for  actual  experi¬ 
mental  predictions,  complicates  a  systematic  comparison  of  numerical  results 
to  theoretical  results.  This  problem  arises  because  current  theoretical  formu¬ 
lations  demand  a  homogeneous  environment  while  the  barium  cloud  model  con¬ 
sists  of  a  striated  cloud-like  formation  with  free  space  to  either  side — an 
inhomogeneous  environment.  For  this  reason,  a  number  of  different  simulations 
have  been  utilized,  with  and  without  the  mean  barium  cloud,  and  are  presented 
in  the  succeeding  sections. 
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4.2 


SIGNAL  PROPAGATION 


To  model  the  proposed  geometry  of  the  PLACES  experiment,  a  single 
phase  screen  representation  of  the  barium  cloud  is  used.  A  plane  v\?ave  is 
then  propagated  through  the  phase  screen  to  a  receiver  plane  100  km  away. 
For  a  thin  scattering  medium  along  the  direction  of  propagation,  this  plane 
wave  geometry  is  equivalent  to  a  spherical  wave  geometry  where  the  random 
medium  lies  between  the  transmitter  and  receiver  such  that 


1  _  JL  +  _L 

100  km  ^  z 

t  r 


(4-4) 


whore  z^  and  z  are  the  distances  of  the  transmitter  and  receiver  from 
the  barium  cloud.  Equation  4-4  is  simply  the  result  of  the  added  angular 
spreading  of  the  spherical  wave  (Reference  4-2). 


4.2.1  Propagation  Results 


Detailed  examples  of  the  numerical  propagation  results  arc  presented 
in  this  subsection.  For  Esther  at  150  Mllz,  the  signal  amplitude  and  phase 
are  shown  as  a  function  of  propagation  distance  from  the  barium  cloud.  The 
effect  of  the  mean  barium  cloud  is  shown  as  well  as  the  effect  of  inner  scale 
size  variation.  Results  for  the  received  signal  amplitude  at  100  km  propaga¬ 
tion  distance  arc  presented  for  both  Esther  and  Fern  at  frequencies  of  150, 
200,  and  250  MHz. 


I'igure  4-1  shows  a  phase  screen  representation  of  Esther  where 
the  mean  barium  cloud  part  of  the  model  has  been  deliberately  omitted.  This 
representation  was  generated  using  the  PSD  parameters  listed  in  Table  4-1 
along  with  an  assumed  inner  scale  of  10  m.  The  phase  shown  has  both  positive 
and  negative  portions  and  represents  a  portion  of  a  zero  mean  random  process. 


*  Current  experimental  plans  call  *'or  transmission  from  a  rocket  borne 
transmitter  through  a  striated  barium  cloud  to  a  ground  based  receiver. 
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Figure  4--1.  Phase  screen  representation  of  Esther  at  150  MHz 
with  no  deterministic  barium  cloud. 
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Note  that  the  MPS  propagation  simulation  is  two-dimensional  with 
no  variation  perpendicular  to  the  page.  This  is  often  a  reasonable  assumption 
in  ionospheric  propagation  cases  where  there  is  little  variation  parallel  to 
the  magnetic  field  direction.  A  wave  (initially  plane  as  it  enters  the 
striated  region)  is  then  propagated  numerically  downwards  to  the  receiver  by 
the  use  of  the  Fresnel-Kirchhoff  integral  equation.  As  the  wave  propagates, 
it  is  possible  to  trace  the  signal  development. 

Figure  4-2  shows  the  amplitude  (in  decibels  relative  to  the  mean 
amplitude)  and  phase  resulting  after  propagating  a  distance  of  1  km  from  the 
striation  location.  At  this  point  severe  amplitude  fluctuations  have  already 
developed.  At  the  1  km  propagation  distance,  however,  the  signal  phase  still 
resembles  the  phase  at  the  striation  location,  although  there  is  a  tendency 
for  the  dominance  of  positive  values  of  phase.  These  increasing  signal  phase 
values  correspond  to  increasing  signal  time  delay  due  to  propagation  along 
longer  paths  resulting  from  angular  scattering  by  the  random  striations. 

Figures  4-3  to  4-7  show  the  development  of  the  signal  amplitude 
and  phase  as  the  pi'opagation  distance  gradually  increases  to  100  km.  In  the 
phase  plots,  the  values  shown  correspond  to  additional  phase  increments  after 
the  free  space  phase  delay  (e  o  )  is  removed.  Thus  only  the  effects  of 
the  striations  are  shown.  When  the  propagation  distance  is  small,  energy 
reaching  a  particular  point  can  originate  only  from  a  relatively  small  portion 
of  the  striations  and  hence  cause  only  a  relatively  small  time  delay.  As 
the  propagation  distances  increases,  energy  from  a  larger  portion  of  the 
barium  cloud  can  reach  a  particular  observer  by  angular  scattering.  Larger 
time  delay  and  the  associated  greater  phase  advance  shown  in  the  figures 
arise  because  of  these  angular  scattering  effects. 

Also  note  the  increased  spreading  of  the  signal  due  to  angular 
scattering  which  becomes  more  important  as  the  propagation  distance 
increases.  This  is  particularly  evident  in  the  amplitude  behavior  at 
the  larger  propagation  distances  shown. 
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Figure  4-2.  Signal  aniplitude  and  phase  at  a  propagation  distance  of  1  km 
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Figure  4-3.  Signal  amplitude  and  phase  at  a  propagation  distance  of  4  km 
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It  is  seen  in  I’igures  4-2  to  4-7  that  the  signal  phase  starts 
at  zero  radians  at  an  abscissa  value  of  zero  on  the  plots  but  does  not  return 
to  zero  at  a  distance  of  30  km.  This  behavior  is  caused  by  rapid  phase 
changes  during  deep  fades  which  have  faster  variation  than  allowed  for  in  the 
grid  sampling.  Wlicn  the  phase  is  reconstructed  from  the  quadrature  components 
in  the  simulation,  there  is  a  finite  probability  of  missing  27T  changes  in  the 
phase  and  these  2it  shifts  are  reflected  by  deviations  from  zero  radians  at  the 
end  of  the  grid.  Since  those  rapid  pliase  changes  occur  during  deep  fades 
where  tliere  is  relatively  little  signal  power,  tliey  have  no  significant 
effect  on  any  results  presented  here. 

I'igure  4-8  shows  the  signal  amplitude  and  pliase  at  a  propagation 
distance  of  100  km  for  a  deterministic  or  mean  barium  cloud  only.  As  can 
be  seen,  the  deterministic  model  of  Usther  (a  gaussian  with  no  striations) 
causes  only  a  slight  defocusing  even  at  a  frequency  of  150  Mllz. 

IMgures  4-0  to  4-14  show  tiie  received  signal  amplitude  at  a  propa¬ 
gation  distance  of  100  km.  These  six  figures  utilize  the  full  barium  cloud 
models  of  I'.sthcr  and  I'ern  ns  tabulated  in  Talile  4-1.  That  is,  both  mean 
cloud  and  striations  are  used  to  obtain  those  results.  For  these  figures 
no  explicit  inner  scale  was  used  in  the  MPS  phase  PSO,  so  that  tlie  effective 
inner  scale  size  was  0.58  m  as  determined  from  the  grid  point  spacing  (a 
30  km  grid  with  10384  points).  A  comparison  of  Figure  4-7  (inner  scale  of 
10  111)  and  I'igure  t-9  (inner  scale  of  0.58  m)  shows  that  the  effect  of  de¬ 
creasing  the  inner  scale  size  is  to  increase  the  amount  of  angular  scattering. 

Figures  4-0  to  4  11  show  the  recced  signal  amplitude  for  Fsther 

at  frequencies  of  150,  200  and  250  Mllz.  The  next  three  figures  show  received 

signal  amplitude  for  Fern  for  the  same  three  frequencies.  It  is  apparent 

from  tlie  figures  that  Fsther  poses  a  more  sex  ere  propagation  environment  than 

docs  Fern.  Of  course  the  figures  simj)!)'  reflect  the  barium  cloud  mxxleling 

parameters  listed  In  Tai'lc  4-1  where  o.  for  Fsther  is  almost  twice  as 
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large  as  for  Fern.  From  the  figures,  it  is  evident  that  the  large 

gaussian  shape  used  to  model  the  deterministic  portions  of  the  barium  cloud 
causes  a  large  defocusing  region  below  the  center  of  the  cloud.  This  de- 
focusing  is  stronger  for  Fern  than  for  Esther  because  the  Fern  cloud  is 
narrower  with  a  larger  mean  phase  variation.  That  is,  from  Table  4-1,  r^ 
is  smaller  and  (j)^  is  larger  for  Fern  than  for  Esther. 

Work  is  currently  in  progress  to  relate  the  propagation  results 
described  above  to  analytical  formulations  for  mean  time  delay,  time  delay 
jitter,  and  signal  temporal  broadening.  These  results  will  be  presented  in 
a  subsequent  report. 

4.2.2  Frequency-Selective  Propagation  Results 

In  this  subsection,  some  examples  of  the  PN  code  correlation 
function  development  with  time  are  presented.  The  effect  of  inner  scale 
variation  is  shown  pictorially  as  well  as  the  effect  of  changing  the  carrier 
frequency  from  150  to  200  ^0^z. 

As  illustrated  in  Section  5  of  this  report,  the  PN  code  correlation 
function  provides  a  measure  of  the  impact  of  frequency-selective  fading  on 
the  operation  of  the  X-set  code  tracking  loop.  In  the  absence  of  frequency- 
selective  effects,  the  envelope  of  the  PN  code  correlator  output  is  triangular 
in  shape  and  the  code  loop  has  little  difficulty  tracking  the  peak,  provided 
sufficient  signal  power  is  available. 

Figure  4-15  shows  a  three-dimensional  plot  of  the  envelope  of  the 
PN  code  correlator  output  for  Esther  at  48  minutes  after  release  at  a  carrier 
frequency  of  150  MHz.  The  PN  code  chip  rate  (10.23  Mbps)  is  much  larger  than 
the  fading  rate  here  so  that  propagation  conditions  remain  unchanged  during  the 
many  chip  durations  associated  with  the  calculation  of  any  one  of  the  PN  code 
correlator  output  traces  shown.  Subsequent  calculations  of  the  correlator 
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Figure  4-15.  PN 
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output,  corresponding  to  successive  points  in  the  MPS  grid,  are  shown  plotted 
behind  previous  calculations  in  the  figure.  Figure  4-15  shows  the  received 
PN  code  correlation  envelope  as  a  function  of  time  delay  in  code  chips  (1 
chip  duration  is  97.75  nsec).  The  plot  shows  101  successive  correlator 
output  traces  corresponding  to  101  successive  points  in  the  MPS  grid  sepa¬ 
rated  by  1.831  m.  At  an  aircraft  velocity  of  250  m/sec  these  101  points 
would  correspond  to  a  time  interval  of  about  0.7  seconds.  The  time  domain 
signal  was  reconstructed  from  16  frequencies  propagated  over  a  bandwidth  of 
20.46  NOIz  so  that  Af  is  1.27875  MHz.  The  unambiguous  period  of  the  recon¬ 
structed  output  waveforms  is  then  1/Af=  7.82x10  ^  sec  (8  chips).  The 
complete  barium  cloud  model  listed  in  Table  4-1  for  Esther  at  48  minutes 
was  used  to  obtain  the  results  shown  in  Figure  4-15.  The  effective  inner 
scale  here  is  0.58  m.  The  101  successive  points  shown  in  the  figure  are 
located  almost  directly  beneath  the  center  of  the  barium  cloud,  i.e.,  near 
the  15  km  distance  in  Figui*e  4-9.  (Figure  4-9  shows  one  of  the  16  frequencies 
used  in  the  OFT  algorithm  to  reconstruct  the  time  domain  signal.)  The  mean 
total  elccti’on  content  (TEC)  through  the  center  of  this  cloud  causes  a  time 
delay  of  about  2  chips  at  150  MHz.  The  delay  due  to  mean  TEC  is  further 
increased  by  frequency- selective  scattering  as  shouTi  in  Figure  4-15. 
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Figure  4-16  shows  the  same  portion  of  the  PN  code  correlator 
output  envelope  as  the  previous  figure,  except  the  inner  scale  is  changed 
to  10  m.  .‘\gain  the  effect  of  the  larger  inner  scale  size  is  to  decrease 
the  amount  of  angular  scattering  and  consequently  produce  smoother  time 
domain  waveforms.  As  discussed  in  Section  3,  the  frequency-selective  propa¬ 
gation  envii’onment  acts  to  distort,  lengthen  and  delay  the  signal  as  shown 
in  the  figure. 

In  both  Figures  4-15  and  4-16  the  abscissa  scales  include  the 
additional  signal  delay  due  to  propagation  through  the  mean  barium  cloud 
as  noted  above.  The  selection  of  the  u.iambiguous  plotting  interval  (DFT 
window)  ^-'or  these  two  figures  is  somewhat  erroneous.  The  signal  energy 
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appearing  at  the  left  hand  side  of  each  plot  should  have  been  plotted  at  the 
right  hand  side  so  that  all  energy  is  delayed  relative  to  the  onset  of  the 
PN  code  correlator  output. 

Figures  4-17  and  4-18  show  a  comparison  of  the  PN  code  correlator 
outputs  for  carrier  frequencies  of  150  and  200  MHz  for  Esther  at  48  minutes 
after  release.  In  both  of  these  cases  the  propagation  simulation  included 
only  the  effects  of  the  striations,  with  no  deterministic  barium  cloud.  Thus 
the  abscissa  scales  sho\>m  in  these  two  figures  do  not  include  the  mean  TEC 
delay.  A  value  of  10  m  was  used  for  the  inner  scale.  Again  the  101  MPS 
points  shown  correspond  to  the  same  locations  approximately  beneath  the  center 
of  the  barium  cloud  striations.  (See  Figure  4-1  for  the  phase  screen  repre¬ 
sentation  of  the  barium  cloud  striations  at  150  MHz.)  For  the  two  figures, 

32  frequencies  were  used  over  a  bandwidth  of  20,46  MHz  with  a  smaller  Af 
and  a  larger  unambiguous  time  interval  of  16  chips.  This  finite  bandwidth, 
which  corresponds  to  the  first  null-to-first  null  width  of  the  PN  signal 
spectrum,  causes  slight  rounding  in  the  triangular  shape  of  the  autocor¬ 
relation  function  and  a  power  reduction  of  0.89  dB  at  the  peak.  A  very  small 
amount  of  additional  rpreading  is  introduced  by  mean  dispersive  effect;  due 
to  an  assumed  total  electron  content  of  5  x 10^  electrons/m  .  Both  effects 
are  included  in  the  two  figures,  although  lIiO  mean  TEC  time  delay  is  not 
included. 

The  effect  of  increasing  the  carrier  frequency  and  thus  reducing 
the  angular  scattering  is  evident  from  tlie  figures.  At  .i0  MHz  energy  is 
delayed  up  to  about  5  chips  as  shown  in  Figure  4-17;  at  a  carrier  frequency 
of  200  MHz  very  little  energy  is  delayed  more  than  3  to  4  chips. 
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Figure  4-17.  pN  code  correlator  output  envelope  for  Esther  at  -i  carrier 
frequency  of  150  MHz. 


Figure  4-'18.  PN  code  correlator  output  envelope  for  Esther  at  a  carrier 
frequency  of  200  MHz. 


4.3 


X-SET  PERFORMANCE  RESULTS 


Results  presented  in  Reference  4-3  illustrate  the  performance  of 
the  X-set  receiver  under  both  normal  and  disturbed  nonselective  signal  con¬ 
ditions.  In  this  section  the  X-set  simulation,  modified  as  described  in 
Section  3  to  include  frequency-selective  effects,  is  exercised  with  MPS 
calculations  of  signal  structures  to  investigate  receiver  performance  under 
disturbed,  frequency-selective  propagation  conditions.  A  variety  of  X-set 
simulations  have  been  performed  for  both  Esther  and  Tern  at  frequencies  from 
150  to  250  MHz.  The  examples  shown  here  are  confined  to  the  32  frequency 
MPS  calculations  of  Esther  at  150  and  200  MHz,  portions  of  which  are  shorn 
in  Figures  4-17  and  4-18.  No  deterministic  barium  cloud  is  present  in  these 
examples,  the  inner  scale  is  10  m  and  the  striations  occupy  only  the  center 
portion  of  the  MPS  grid  (as  shown  in  Figure  4-1). 

Of  the  entire  30  km  long  MPS  grid,  only  a  small  portion  (3.75  km) 
was  available  for  utilization  in  the  X-set  receiver  simulation.  This  3.75 
km  length  extended  from  11.25  km  to  15  km  on  the  abscissa  scale  shown  in 
Figure  4-1  and  consequently  included  the  region  near  the  center  of  the 
striated  cloud  where  angular  scattering  effects  are  most  severe.  For  use 
in  the  receiver  simulation,  a  velocity  of  250  m/sec  is  assumed  and  the  3.75 
km  stretch  of  data  is  converted  to  a  15  second  time  period.  This  period  is 
too  short  to  accumulate  statistics  on  receiver  performance.  However,  the 
variety  of  shoi’t  simulation  runs  that  have  beei.  performed  do  provide  valuable 
insight  into  the  manner  in  which  frequency-s.,1 active  propagation  disturbances 
impact  X-set  receiver  operation. 

4.3.1  Receiver  Design  Parameters 

Except  where  otherwise  noted,  results  presented  in  this  section 
are  obtained  using  receiver  design  parameter  values  listed  in  Table  4-2. 

These  parameters  are  discussed  in  Reference  4-3  and  many  results  are  presented 
there  using  this  baseline  set  of  design  values. 


96 


I* 

\ 


t-' 

( 


S'. 


As  shown  in  Table  4-2,  the  third-order  Costas  loop  bandwidth  is 
21.3  Hz,  the  second-order  AFC  loop  bandwidth  is  3.3  Hz  and  the  first-order 
code  loop  bandwidth  is  0.1  Hz.  Doppler  (velocity)  aiding  from  the  carrier 
tracking  loop  is  applied  to  the  code  tracking  loop.  These  and  the  other 
design  values  listed  in  the  table  are  all  inputs  to  our  X-set  simulation 
program. 

With  one  important  exception,  the  values  shown  in  Table  4-2  are 
obtained  primarily  from  GPS  Phase  I  specification  documents  (References  4-4, 
4-5,  4-6)  and  from  the  Draper  Laboratory  report  by  Stonestreet  (Reference  4-7) 
The  major  exception  is  that  we  have  gone  back  to  the  250/sec  AFC/Costas 
loop  iteration  rate  used  by  Cahn  (Reference  4-8)  for  the  baseline  design, 
instead  of  the  50/sec  rate  stated  in  Reference  4-7.  This  is  done  for  two 
reasons:  (1)  at  the  larger  Costas  loop  bandwidths  of  interest,  the  250/sec 
loop  iteration  rate  provides  better  X-set  performance,  and  (2)  discussions 
with  GPS  JPO  personnel  indicate  that  this  iteration  rate  may  be  I’einstated 
in  Phase  II  or  Phase  III  (Reference  4-9). 

Table  4-2.  Nominal  design  values  used  in  X-set  simulation. 


AFC/COST«S/PN  CPS  X-SET  RECEIVER  OPERrtTIOM  IM  SlIlULftTED  BORIUII  CLOUD  ENUIROMMEMT 

RECEIVER  DESIGN  PARRHETERS 
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UPDATE  RATE 
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PN  CODE  CHIP  RATE  (BPS) 

1 .oasooEre? 

AGC 

TIME  CONSTANT 

(SEC) 

I .eoeoeEiOO 

NAV  DATA  BIT  RATE  (BPS) 

s.eoeooEiOi 
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a 
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- 
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8 

CODE 

UPDATE  TIME 
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5 
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1 
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) 
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CODE 
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AFC 
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3 

AFC 
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a 
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3.300OOEi00 
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AFC 

DAMPING  FACTOR 
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COSTAS  LOCK  TINE  (SEC) 

a.40O00E-0I 
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!OCK  TIME 

(SEC) 

9.60000E-O1 

COSTAS  LOCK  THRESHOLD 
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AFC 
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The  carrier  loop  iteration  rate  is  one  of  the  design  excursions 
investigated  in  several  disturbed  propagation  cases  in  this  section.  Ox  the 
various  AFC/Costas  loop  filter  implementations  tested  previously  at  SO/ sec 
iteration  rate  with  nominal  21.3  Hz  Costas  loop  bandwith,  one  implementation 
stood  out  as  by  far  the  best  choice.  This  implementation  employs  the  "mean- 
value"  algorithm,  which  simply  represents  a  specific  method  of  choosing  loop 
filter  coefficients,  in  conjunction  with  lumping  the  foreground  portion  of 
the  filter  output  and  constraining  the  foreground  delay  to  be  in  t.hx  0  to 
4  ms  range.  With  this  specific  implementation  of  the  Costas  loop,  fairly 
good  X-set  performance  is  obtained  at  50/sec  loop  iteration.  There  is  still 
noticeable  degradation  of  carrier  and  code  tracking  performance  compared  to 
that  obtained  with  250/sec  carrier  loop  iteration,  but  the  degradation  is 
not  nearly  as  serioius  as  with  other  implementations  tested.  Reference  4-3 
contains  detailed  information  on  this  subject. 

The  values  shown  in  Table  4-2  for  the  Costas  and  AFC  lock  detection 
filter  integ-^ation  times  and  threshold  settings  are  somewhat  speculative  on 
our  part.  These  parameters  are  discussed  in  Reference  4-3  and  are  based  in 
part  on  information  in  References  4-7  and  4-8.  The  Costas  lock  detection 
filter  is  used  to  enable  the  AFC  loop  when  Costas  lock  is  lost  and  to  disable 
the  AFC  when  Costas  lock  is  detected.  This  is  the  automatic  AFC  switching 
mode  assumed  here  for  the  baseline  design  with  250/sec  AFC/ Costas  loop 
iteration  rate.  The  Costas  lock  detection  filter  state  also  determines  the 
AGC  measurement  mode.  The  AGC  operates  coherently  when  Costas  lock  is 
detected  and  noncoherently  when  loss  of  Costas  lock  is  sensed. 

The  AFC  lock  detection  filter  provides  simulation  output  indicating 
loss  or  reacquisition  of  frequency  lock,  but  has  no  internal  effect  on  simu¬ 
lation  operation  here.  In  the  actual  X-set  receiver,  it  is  our  understanding 
that  indication  of  AFC  out-of-lock  might  initiate  a  code  reacquisition  pro¬ 
cedure.  If  this  is  the  case,  the  AFC  lock  indicator  is  an  important  opera¬ 
tional  consideration.  The  AFC  lock  integration  time  shown  in  Table  4-2 
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actually  corresponds  to  a  1.2  sec  interval  because  one  of  every  five  A/D 
samples  straddles  a  data  bit  edge  and  is  not  used  in  the  AFC  lock  detection 
accumulation.  The  1.2  sec  integration  period  is  based  on  results  obtained 
in  previous  investigations  where  it  was  found  that  a  shorter  integration  time 
caused  the  filter  to  indicate  loss  of  lock  during  brief  fades  when  AFC  lock 
was  actually  maintained.  The  integration  interval  should  not  be  too  large 
because  long,  deep  signal  fades  can  indeed  cause  loss  of  AFC  and  code  lock. 

The  AGC  design  value  of  C/N^  shown  in  Table  4-2  is  used  in  the 
simulation  to  scale  input  signal  amplitude  and  noise  samples  to  the  nominal 
design  level.  In  effect,  this  parameter  initializes  the  AGC  gain  adjustment 
in  the  simulation.  To  avoid  possibly  large  AGC  transients  at  the  beginning 
of  each  simulation  trial  here,  the  AGC  is  initialized  to  the  mean  input 
signal  level  at  the  start  of  each  case.  Thereafter,  the  AGC  loop  operates 
in  a  normal  manner. 

A  detailed  description  of  the  GPS  navigation  data  formatting  ;nd 
error-detection  encoding  is  contained  in  Reference  4-6.  The  data  bit  7..' 
is  50  bps.  complete  data  message  (called  a  frame)  consists  of  150C  )  's, 
and  is  divided  into  five  subframes.  Each  subframe  consists  of  ten  .50 -bit 
words.  Thus  a  word  is  transmitted  in  0.6  sec,  a  subframe  in  6  sec,  and  a 
complete  data  frame  in  30  sec.  Data  tra>‘smission  is  continuous,  rnci  frames 
are  therefore  repeated  at  30  sec  intervals. 

All  data  is  encoded  using  a  (32,26)  Hamming  code,  wh’ch  is  a  single 
error-correction,  double  error-detection  code.  The  X-set  receiver  simulation 
incorporates  the  GPS  data  formatting  and  encoding  described  in  Reference  4-6. 
The  error- correct  ion  capability  of  the  Hamming  code  is  very  limited  and  is 
not  used.  However,  the  error-detectirn  capability  is  quite  useful.  Also, 
the  (32,26)  Hamming  code  provides  resolution  of  the  tt  phase  ambiguity 
inherent  in  Costas  phase  tracking. 
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4.3.2 


X-Set  Performance  at  150  MHz 


A  few  graphical  examples  of  X-set  receiver  simulation  results  are 
presented  to  provide  insight  into  the  effects  of  •'ignal  scintillation  on 
receiver  operation. 

Figure  4-19  shows  the  input  signal  amplitude  and  phase  over  a 
15  second  simulation  time  interval  at  a  carrier  frequency  of  150  NiHz.  Signal 
amplitude  is  plotted  in  decibels  relative  to  the  mean  signal  level.  The  cor¬ 
responding  signal  phase  is  plotted  in  radians..  The  time  scale  here  results 
from  a  250  m/sec  conversion  of  the  distance  scale  used  in  the  numerical  propa¬ 
gation  algorithm.  This  signal  realization  is  the  sixteenth  member  of  the 
thirty-two  frequency  comb  used  as  receiver  input  to  constr>'.ct  the  received 
time  domain  signal  at  the  PN  code  correlator  output. 

Figures  4-20  through  4-22  compare  the  resulting  phase,  frequency, 
and  code  tracking  performance  of  the  simulated  X-set  receiver  for  two 
different  values  of  mean  carrier  power-to-noise  density  ratio  (C/N^) ,  35 
and  32  dB-Ilz. 

Figure  4-20  shows  the  phase  tracking  error  time  history  for  the 
two  values  of  mean  C/N^  for  the  15  second  simulation  period.  This  phase 
tracking  error  is  obtained  in  the  simulation  by  differencing  the  incoming 
value  of  the  signal  phase  with  the  estimated  phase  developed  in  the  AFC/Costas 
loop  and,  as  such,  is  not  available  to  the  receiver  itself.  Three  adjacent 
equilibrium  phase  tracking  states  are  shown  on  the  plot  at  0  and  ±ir  radians. 
The  upper  figure,  with  a  mean  value  of  35  dB-Hz  exhibits  much  fewer 

phase  slips,  and  only  two  brief  periods  of  loss  of  phase-lock  at  about  7.5 
seconds  and  at  about  12  seconds.  In  both  cases,  rhis  temporary  loss  of 
phase-lock  is  associated  witli  simultaneous  large  receiver  estimates  of 
Doppler  frequency  sho^vn  in  the  top  of  Figure  4-21.  The  lower  portion  of 
Figure  4-20,  at  a  mean  C/N^  value  of  32  dB-tlz  is  characterized  by  very 
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rough  phase  tracking  and  two  extended  periods  of  loss  of  phase- lock  which, 
are  caused  by  a  quite  complicated  interaction  of  the  AFC/Costas  loop  and 
the  code  tracking  loop  discussed  below.  Overall,  the  phase  tracking  per¬ 
formance  for  this  highly  disturbed  barium  environment  can  be  characterized 
as  quite  noisy  with  many  phase  slips.  This  poor  phase  tracking  is  reflected 
in  poor  data  demodulation  performance  and  hence  large  word  error  rates. 

Figure  4-21  shows  a  comparison  of  the  estimated  Doppler  “frequency 
developed  in  the  X-set  AFC/Costas  loop  filter.  The  lower  figure,  for  a  mean 
C/N^  value  of  32  dB-Hz  exhibits  much  greater  frequency  deviations  (note  that 
the  ordinates  are  different  for  the  upper  and  lower  plots)  and  rougher  track. 
In  particular,  note  the  large,  negative  frequency  deviation  occurring  between 
5  and  6  seconds  simulation  time.  This  large  excursion  indicates  that  the 
AFC 'Costas  loop  is  out  of  phase  lock  (as  shown  in  the  previous  figure)  and 
also  temporarily  out  of  frequency  lock.  This  in  turn  influences  the  code 
tracking  loop,  which  is  Doppler  aided  from  the  AFC/Costas  loop. 

Figure  4-22  shows  a  comparison  of  the  code  tracking  performance 
of  the  X-set  at  the  two  values  of  mean  C/N^.  Again  note  that  the  ordinate 
scales  are  different.  As  will  be  shown,  the  X-set  at  35  dB-Ilz  C/N^  ade¬ 
quately  tracks  the  peak  of  the  PN  code  correlation  envelope.  At  about  12 
seconds  in  the  upper  plot  and  around  10  seconds  in  the  lower  plot  there  is 
a  rapid  change  in  the  estimated  code  delay  caused  by  Doppler  aiding  from  the 
AFC/Costas  loop.  Note  the  corresponding  large  frequency  excursions  at  these 
times  in  Figure  4-21.  On  the  bottom  plot  of  Figure  4-22,  similar  Doppler 
"aiding"  to  the  code  loop,  at  about  5  seconds  simulation  time,  acts  to  drive 
the  code  delay  estimate  off  the  main  correlation  peak  to  an  extended  sidelobe. 
The  code  loop  then  ti’acks  this  sidelobe  for  about  5  seconds  before  reacquiring 
the  main  peak  at  about  10  seconds  simulation  time. 
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Figure  4-23  shows  a  comparison  of  the  receiver  code  correlator 
amplitude  for  the  two  cases.  The  amplitude  scale  is  decibels  relative  to 
the  nominal  correlator  output  level.  (Under  ideal  conditions,  and  perfect 
code  tracking,  the  code  correlator  amplitude  takes  on  a  value  of  0  dB  cor¬ 
responding  to  the  peak  of  the  PN  code  autocorrelation  function  envelope.) 
Except  for  the  period  between  5  and  10  seconds,  the  two  plots  are  almost 
identical.  As  discussed  above,  the  differences  are  caused  by  tracking  of  a 
sidelobe  by  the  code  loop  during  this  interval.  During  the  sidelobe  tracking 
period,  less  power  appears  at  the  correlator  output  and  consequently  the 
X-set  performs  as  if  the  mean  C/N^  were  reduced. 

Figures  4-24  through  4-26  show  successive  plots  of  the  PN  code 
correlator  output  envelope  for  three  brief  periods  in  the  15  second 
simulation.  Each  plot  shows  about  0.73  seconds  of  simulation  time.  Figure 
4-24  shows  the  period  between  6.60  sec  and  7.33  sec.  Figure  4-25  shows  the 
period  from  10.26  sec  to  10.99  sec  and  Figure  4-26  shows  the  interval  from 
13.93  sec  to  14.66  sec.  Note  the  increasingly  large  distortion  of  the  PN 
code  correlation  envelope  as  the  simulation  time  increases  to  15  seconds 
(Figure  4-26),  at  which  point  the  observer  is  almost  directly  below  the 
center  of  the  barium  cloud  striations  and  thus  subject  to  angular  scattering 
arising  from  striations  on  either  side  of  the  cloud  center.  As  the  distance 
from  the  point  below  the  cloud  center  increases  (with  decreasing  simulation 
time)  the  angular  scattering  effects  are  less  severe,  as  may  be  seen  by 
comparing  Figures  4-24  and  4-26. 

Also  shown  in  Figures  4-24  through  4-26  as  a  series  of  arrows  is 
the  X-set  estimated  time  delay  as  it  is  tracked  in  t^ e  receiver  delay-lock 
loop  for  the  two  values  of  C/N^  of  35  and  32  dB-Hz.  These  estimated  time 
delays  in  the  code  loop  are  the  same  delay  estimates  plotted  in  Figure  4-22. 
In  Figure  4-24,  for  the  35  dB-Hz  C/N^  case  the  code  loop  can  be  seen  to  track 
the  first  major  peak  of  the  PN  code  envelope  quite  well  at  a  time  delay  of 
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Esther  at  150  MHz:  PN  code  envelope  and  code  tracking  history 
for  simulation  time  from  6.60  to  7.33  seconds. 


Esther  at  150  MHz:  PN  code  envelope  and  code  tracking  history 
for  simulation  time  from  10.26  to  10.99  seconds. 


-26.  Esther  at  150  MHz:  PN  code  envelope  and  code  tracking  history 
for  simulation  time  from  13.93  to  14.66  seconds. 


about  0.3  chips  throughout  the  0.73  seconds  of  simulation  time  shown.  For 
the  32  dB-Hz  C/N^  case,  erroneous  Doppler  aiding  from  the  AFC/Costas  loop 
has  caused  the  code  loop  to  lock  onto  a  sidelobe  at  around  1.3  to  1.5  chips 
time  delay. 

In  Figure  4-25,  at  a  simulation  time  interval  from  10.26  to  10.99 
sec,  the  X-set  operating  at  a  mean  C/N^  value  of  35  dB-Hz  is  tracking  the 
peak  quite  well.  At  C/N^  of  32  dB-Hz,  the  receiver  is  still  recovering  from 
the  extended  5  second  period  during  which  the  code  loop  was  tracking  a  cor¬ 
relation  sidelobe. 

A  little  later  in  the  15  second  simulation,  the  code  loop  reacquires 
the  peak  and  each  of  the  two  values  of  C/N^  give  the  same  code  tracking  per¬ 
formance.  During  the  brief  interval  sliown  in  Figure  4-26,  the  X-set  operation 
at  32  dB-Hz  C/N^  is  comparable  to  the  operation  at  35  dB-Hz,  at  least  as  far 
as  frequency  and  code  tracking  is  concerned. 

Other  measures  of  X-set  performance  are  also  degraded  at  the  lower 

value  of  mean  C/N^.  During  the  15  second  simulation  interval,  25  30-bit 

words  are  received  and  decoded  using  the  (32,26)  Hamming  code  described  in 

Reference  4-6.  For  a  C/N  value  of  35  dB-Hz,  16  of  the  25  words  are  received 

0 

with  one  or  more  bit  errors,  giving  a  word  error  rate  of  0.64.  For  the 
lower  C/N^  value  of  32  dB-llz,  only  one  word  is  received  free  of  errors,  for 
a  word  error  rate  of  0.96.  For  35  dB-Hz  mean  C/N^  the  rms  time-delay  error 
is  J.37  chips  corresponding  to  a  position  error  of  10.8  meters  (1  ciiip  of 
the  10.23  Mbps  P  code  equals  approximately  29.305  meters);  for  the  32  dB-Hz 
mean  C/N^  case  the  rms  time-delay  error  is  0.92  chips  corresponding  to  a 
position  error  of  26.96  meters. 

Similar  15  second  simulations  were  performed  with  the  50/sec 
AFC/Costas  loop  iteration  rate  for  mean  C/N^  values  of  44,  38,  and  32  dB-Hz. 
For  this  lower  iteration  rate  with  worse  phase  and  frequency  tracking  per¬ 
formance  than  the  250/sec  rate,  loss  of  frequency  lock  in  the  AFC/Costas 
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loop  was  found  to  drive  the  code  loop  out  of  lock  because  of  the  Doppler 
aiding  applied  in  the  code  loop.  All  three  values  of  C/N^  exhibited 
this  same  loss  of  frequency  and  code  lock  condition.  No  effort  was  made  to 
investigate  code  reacquisition  in  this  work. 

4.3.3  X-Se^  Performance  at  200  MHz 

Three  15  second  simulations  at  a  carrier  frequency  of  200  MHz  were 
performed  for  the  same  barium  model  of  Esther  (without  the  deterministic 
cloud)  as  those  described  in  Section  4.3.2.  X-set  performance  was  investi¬ 
gated  for  a  value  of  mean  C/N^  of  35  dB-Hz  with  AFC/Costas  iteration  rates  of 
250/sec  and  50/sec,  and  for  a  mean  C/N^  value  of  32  dB-Hz  with  an  iteration 
rate  of  250/sec.  In  general,  code  tracking  performance  was  quite  good,  with 
no  problems  for  the  higher  iteration  rate  at  either  value  of  C/N^.  For  the 
slower  50/sec  AFC/Costas  iteration  I'ate,  incorrect  Doppler  aiding  from 
the  frequency  tracking  loop  drove  the  code  loop  briefly  out  of  lock  twice, 
but  recovery  occurred  rapidly  each  time.  For  the  250/sec  iteration  rate 

the  word  error  rates  were  0.48  and  0.64  for  C/N  of  35  dB-Iiz  and  32  dB-ilz, 

o 

respectively.  For  the  50/sec  iteration  rate  at  a  mean  C/N^  value  of  35  dB-Hz, 
the  word  error  rate  was  0.68. 

4.3.4  Summary  of  X-Set  Performance 

In  addition  to  the  results  discussed  in  the  preceding  subsections, 
simulations  also  have  been  performed  for  carrier  frequencies  of  150,  200, 
and  250  MHz  for  the  full  barium  models  of  both  Esther  and  i orn  (deterministic 
barium  cloud  included)  with  an  effective  inner  scale  of  0.58  m.  For  these 
cases,  the  major  effect  of  the  deterministic  cloud  is  to  cause  large  scale 
defocusing  which  acts  to  reduce  the  signal  level  in  the  central  region 
immediately  below  the  center  of  the  cloud.  (See  Figures  4-8  through  4-14.) 
This  defocusing  is,  of  course,  equivalent  to  a  reduction  in  the  mean  C/N^ 
level  and  causes  a  resulting  degradation  in  X-set  performance.  For  these 
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six  simulations  and  the  previously  described  simulation  results,  it  is  found 
that  loss  of  code  lock  or  degraded  code  tracking  was  usually  caused  by 
erroneous  Doppler  aiding  from  the  AFC/Costas  loop.  Even  with  good  code 
tracking,  Costas  phase  tracking  and  data  demodulation  is  severely  degraded 
in  the  barium  environment  with  word  error  rates  for  the  30-bit  navigation 
data  word  found  to  be  as  lai'ge  as  0.6  to  1.0. 

In  the  operational  Global  Positioning  System,  large  error  rates 
may  be  tolerable  because  of  redundant  transmission  of  data  frames  every  30 
seconds.  Perhaps  a  more  important  consideration  is  the  undetected  error 
rate.  With  the  GPS  (32,26)  Hamming  code  and  receiver  parity  check  algorithm, 
very  few  words  with  errors  pass  the  parity  check.  Hence  the  undetected 
error  rate  is  quite  low.  Of  the  350  words  transmitted  during  210  seconds 
of  simulation  time  here  only  one  word  incorrectly  passed  the  receiver 
parity  check. 
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SECTION  5 

SUMMARY  AND  RECOMMENDATIONS 


In  this  work,  results  have  been  presented  i'nr  the  performance  of 
the  GPS  X-set  receiver  operating  in  a  striated  barium  environment  at  carrier 
frequencies  from  150  to  250  MHz,  These  results  were  obtained  as  a  simulation 
of  an  airborne  receiver  passing  beneath  a  barium  cloud  and  flying  perpendicular 
to  the  projection  of  the  striations  at  a  velocity  component  of  250  m/sec.  The 
barium  cloud  striations  were  represented  by  a  statistical  model  obtained  from 
the  analysis  of  back-propagation  results  from  the  STRESS  barium  release 
experiment. 

The  receiver  simulation  results  presented  here,  together  with 
earlier  studies  of  X-set  performance  in  a  disturbed  environment  (Reference 
5-lj,  give  reliable  predictions  of  X-set  operation  in  both  nonselective  and 
severe  frequency-selective  propagation  environments.  Detailed  simulation 
examples,  showing  the  complex  interaction  between  the  coupled  X-set  AFC/ 

Costas  and  code  tracking  loops,  have  been  presented.  At  a  carrier  frequency 
of  150  I'Qlz,  the  X-set,  with  250/sec  AFC/Costas  iteration  rate,  performs  the 
critical  function  of  tracking  the  PN  code  quite  well.  However,  erroneous 
Doppler  aiding  from  the  frequency  tracking  loop  can  force  loss  of  code  lock 
for  mean  carrier  power-to-noise  density  ratios  below  about  35  dB-Hz.  Word 
error  rates  for  the  30-bit  navigation  data  words  are  found  to  be  as  large 
as  0.6  to  1.0,  even  with  good  code  tracking. 

Because  of  problems  associated  with  the  transmittal  of  large 
amounts  of  propagation  results  required  for  the  frequency-selective  receiver 
simulation,  it  was  difficult  tc  obtain  X-set  simulation  results  of  sufficiently 
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long  duration  for  statistical  significance.  This  difficulty  has  been 
alleviated  and  future  work  should  include  the  compilation  of  more  results 
for  X-set  performance.  Further  predictions  should  be  made  for  a  carrier 
frequency  of  100  MHz  with  consideration  of  the  planned  PLACES  rocket  geometry. 

Work  is  currently  in  progress  to  relate  the  propagation  calcula¬ 
tions  presented  here  to  analytical  results  for  mean  time  delay,  time  delay 
jitter,  and  broadening  of  the  received  waveform.  These  issues  will  be 
addressed  in  a  separate  report. 

REFERENCES 

5-1.  Bogusch,  R.  L.,  and  D.  L.  Knepp,  Propagation  Effects  on  GPS  Receiver 
Model  X.  Mission  Research  Corporation,  January  1978 ’(Unpubl ished) . 


116 


1 


I 

I 

%-  'K 


i  '1-v 


i*gjE>-Mi«yi«j|f>apc|--t-niT - -  '*  -  ‘ _ Jlr...£^Ls _ _ _ _ _  _ <... 


1  ■■■twW^nwwv<rlf^tfn«gn^^ 


DISTRIBUTION  LIST 


Assistant  Secretary  of  Defense 
Comm.,  Cmd,  Cont. ,  &  Intell 

ATTN:  Dir  of  Intelligence  Systems,  J. 
ATTN:  C3IST&CCS,  M.  Epstein 


DEPARTMENT  OF  THE  NAVY  (Continued) 
Na«1  Ocean  Systems  Center 


Babcock 


Assistant  to  the  Secretary  of  Defense 
Atomic  Energy 

ATTN:  Executive  Assistant 

Defense  Advanced  Rsch.  Proj,  Agency 
ATTN:  TIO 

Defense  Communications  Engineer  Center 


ATTN 

ATTN 

ATTN 


Code  R410, 
Code  R720,  0. 
Code  R123 


Craighill 
Worthington 


Defense  Nuclear  Agency 


ATTN 
ATTN 
cy  ATTN 


DDST 

STVL 

RAAE 

TITL 


field  Comnand 
Defense  Nuclear  Agency 
AfTN:  FCPR 

Field  Command 
Defense  Nuclear  Agency 
Livonrare  Division 
AT  IN:  FCPR!. 

U.ndersecrelary  of  Defense  for  Rsch.  &  Engrg. 
ATTN;  Strategic  &  Space  Systems  (OS) 

WWMCCS  System  Engineering  Org. 

ATTN:  R.  Crawford 

DEPARTMENT  OF  THE  ARMY 

BHD  Advanced  Technology  Center 
Department  of  the  Army 


ATTN 

ATTN 

ATTN 


ATC-R.  D. 
ATC-0,  W. 
ATC-T,  H. 


Russ 

Davies 

Capps 


Harry  Diamond  Laboratories 
Department  of  the  Arriy 
ATTN:  DELHD-N  P 

U.S.  Army  Nuclear  &  Chemical  Agency 
ATTN:  Library 

U.S.  Army  Satellite  Comm.  Agency 
ATTN:  Document  Control 

DEPARTMENT  OF  THE  NAVY 

Joint  Cruise  Missile  Project  Office 
Department  of  the  Navy 
ATTN:  JCM-G-70 


ATTN 
ATTN 
3  cy  ATTN 


Code  532,  J.  Bickel 
Code  5322,  M.  Paulson 
Code  5324,  W.  Moler 


Naval  Research  Laboratory 


ATTN 

ATTN 

ATTN 

ATTN 


Code  6700,  T.  Coffey 
Code  7500,  B.  Wald 
Code  6780,  S.  Ossakow 
Code  7550,  J.  Davis 


Naval  Surface  Weapons  Center 
ATTN:  Code  F31 

Office  of  Naval  Research 


ATTN 

ATTN 

ATTN 


Code  465 
Code  421 
Code  420 


Strategic  Systems  Project  Office 


t 

ATTN: 

NSP-43 

\ 

Defense  Technical  Information  Center 

ATTN: 

NSP-2141 

\ 

12  cy  ATTN:  DD 

ATTN: 

NSP-2722,  F.  Wimberly 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Avionics  Laboratory 
ATTN:  AAD,  W.  Hunt 
ATTN:  /iAD,  A.  Johnson 

Air  Force  Geophysics  LabornJory 
ATTN:  bPR-1,  J.  Ult.i,.k 
ATTN:  PHP,  J.  Aarons 
ATTN:  LKB,  K.  Champion 
ATTN:  OPR,  A.  S’air 
ATTN:  OPR,  H.  Gardiner 
ATTN:  PHI,  J.  Buchau 
ATTN:  PHP,  J.  Hulien 

Air  Force  Weapon..  LaboKuory 
ATTN:  Sl'L 
ATTN:  DV: 

Assistant  C.hief  of  Staff 

Studies  &  Analyses 

Departn»nt  of  the  Air  Force 
ATTN:  AF/SASC,  G.  Zank 
ATTN:  AF/SASC,  W.  Adams 

Ballistic  Missile  Office 

Air  Force  Systems  Command 
ATTN:  MNNL 
ATTN:  .MNNH,  M.  Baran 

Deputy  Chief  of  Staff 

Operations  Plans  and  Readiness 

Department  of  the  Air  Force 


ATTN 

ATTN 

ATTN 

ATTN 


AFXOKT 

AFXOXFD 

AFXOKS 

AFXOKCO 


117 


DEPARTMENT  OFItiE  AIR  FORCE  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Deputy  Chief  of  Staff 

BDM  Corp. 

Research,  Development,  S  Acq. 

ATTN: 

L.  Jacobs 

Department  of  the  Air  Force 

ATTN: 

T.  Neighbors 

ATTN:  AFRDS 

ATTN:  AFRDQ 

Berkeley  Research  Associates,  Inc 

ATTN:  AFRDSS 

ATTN: 

J.  Workman 

ATTN:  AFRDSP 

Electronic  Systems  Division 

Comsat  Labs 
ATTN: 

R.  Taur 

Department  of  the  Air  Force 

ATTN: 

G.  Hyde 

ATTN:  DCKC,  0.  Clark 

Electrospace 

Systems,  Inc. 

Electronic  Systems  Division 

ATTN: 

H.  Logs  ton 

Department  of  the  Air  Force 

ATTN:  XRW,  J.  Deas 

ESL,  Inc. 

ATTN: 

J.  Marshall 

Electronic  Systems  Division 

ATTN: 

C.  Prettie 

Department  of  the  Air  Force 

ATTN: 

J.  Roberts 

ATTN:  YSEA 

ATTN:  YSM,  J.  Kobelski 

General  Electric  Company— TEMPO 

Headquarters  Space  Division 

ATTN:  W.  Knapp 
ATTN:  T.  Stevens 

Air  Force  Systems  Command 

ATTN:  W.  McNamara 

ATTN:  E.  Butt 

ATTN:  M.  Stanton 

Strategic  Air  Command 

ATTN:  DAS  1 AC 

ATTN:  0.  Chanaler 

Department  of  the  Air  Force 

ATTN:  DCX 

General  Research  Corp. 

ATTN:  XPFS 

ATTN:  J.  Garbarino 

ATTN:  OOKSN 

ATTN:  J.  Ise,  Jr. 

ATTN:  DCXT 

ATTN:  OCXF 

IBM  Corp. 

ATTN:  NRT 

ATTN:  F.  Ricci 

ATTN:  DCXT,  T.  Jorgensen 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

University  of  Illinois 
ATTN:  K.  Yeh 

Los  Alamos 

Scientific 

Labora tory 

Institute  for  Defense  Analyses 

ATTN:  Document 

Control 

for 

D.  Simons 

ATTN 

J.  Bengston 

ATTN 

H.  Wolfhard 

Sandia  Laboratories 

AHN 

E.  Bauer 

ATTN 

Document 

Control 

for 

Space  Project  Div. 

ATTN 

J.  Aein 

ATTN 

Document 

Control 

for 

Org.  1250,  W.  Brown 

ATTN 

Document 

Control 

for 

3141 

M. 1 .T.  Lincoln  Lab 

ATTN 

Document 

Control 

for 

D.  Oahlgren 

ATTN:  L.  Laughlin 

ATTN 

Document 

Control 

for 

Org.  4241,  T.  Wright 

ATTN:  D.  Towle 

ATTN 

Document 

Control 

for 

D.  Thornbrough 

Mission  Research  Corp. 

EG&G,  Inc. 

ATTN 

F.  Fajen 

ATTN 

Document 

Control 

for 

0.  Wright 

ATTN 

D.  Sowle 

ATTN:  Document 

Control 

for 

J.  Colvin 

ATTN 

R.  Kilb 

ATTN 

R.  Bogusc^' 

DEPARTMENT  OF  DEFENSE 

CONTRACTORS 

ATTN 

R.  Hendrick 

ATTN 

S.  Gutsche 

Aerospace  Corp. 

ATTN 

D.  Sappenfield 

ATTN 

N. 

Stockwell 

ATTN 

D.  Knepp 

ATTN 

V. 

Josephson 

5  cy  ATTN 

Document  Control 

ATTN 

D. 

Olsen 

ATTN 

I. 

Garfunkel 

Mitre  Corp. 

ATTN 

T. 

Salmi 

ATTN 

C.  Callahan 

ATTN 

S. 

Bower 

ATTN 

B.  Adams 

ATTN 

F. 

Morse 

ATTN 

A.  Kynnel 

ATTN 

R. 

Slaughter 

ATTN 

6.  Harding 

Physical  Dynamics,  Inc. 
ATTN:  E.  Fremouw 


R&D  Associates 

ATIN:  L.  Delaney 
ATTN:  B.  Yoon 
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Mitre  Corp 

Rand  Corp. 

ATTN 

M. 

Horrocks 

ATTN:  C.  Crain 

ATTN 

W. 

Foster 

ATTN:  E.  Bedrosian 

ATTN 

W. 

Hall 

ATTN 

J. 

Wheeler 

Science  Applications,  .’nc 
ATTN:  0.  Sachs 

R&D  Associates 

ATTN:  L.  Linson 

ATTN 

R. 

Turco 

ATTN 

C. 

Grei finger 

Science  Applications,  Inc 

ATTN 

'  W. 

Karzas 

AHN:  0.  Divis 

ATTN 

B. 

Gabbard 

ATTN 

W. 

Wright,  Jr. 

SRI  International 

ATTN 

M. 

Gantsweg 

ATTN:  W.  Chesnut 

ATTN 

C. 

MacDonald 

ATTN:  C.  Rino 

ATTN 

F. 

Gilmore 

ATTN;  A.  Burns 

ATTN 

H. 

Ory 

ATTN 

R. 

Lelevier 

> 


